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NIGGLI’S PRINCIPLES OF IGNEOUS PETROGENESIS 
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Received March 25, 1952 


RECENT advances in igneous petrogenesis are based almost entirely on the 
researches made at the Geophysical Laboratory, Washington. The phase 
relationships of various systems that are studied with improved modern 
technique, have given valuable insight into the possible mode of origin of 
rock-forming minerals from silicate melts. Valuable though such knowledge 
is, there remains the fact that the fundamental cause may have to be sought 
for in a sphere that defies experimental analysis. After all experiments are 
directed to determine the physical and chemical conditions that control the 
formation of minerals. But the genesis of minerals and their pattern are 
the natural outcome of a directed dynamic force. 


In the organic world, a plant has its identity determined by the seed, 
and an animal by the egg. It would be difficult to decide how far external 
influences have their effects on the embryo and in what manner. Though 
how a seed gives rise to a particular plant and an egg to a particular animal 
is known, why they should do so is beyond scientific comprehension. A 
mineral too has its pattern determined by forces at work in the melt state. 


Hence our solved problems in the laboratories are not to be considered 
as ends in themselves. There remain undoubtedly problems, the solutions 
of which are to be found only in the realm of speculation. 


Niggli considers that petrogenetic causes are to be traced to reactions 
that take place in the melt. in other words, homogeneous equilibria condi- 
tions. 


The study of petrogenetic and allied problems is usually based on the 
weight percentages of oxides taken as units. According to Niggli, such 
units cannot be anything more than the raw material which require scientific 
and critical treatment to discover hidden values. x grams of MgO do not 
correspond to x grams of Na,O or Al,O,. Therefore “to combine x grams 
of KO and y grams of Na,O to form z grams of alkali oxides, is just as 
meaningless as adding x grams of apples to y grams of potatoes ”’. 


He therefore advocates that atomic and molecular values derived from 
the weight percentages of oxides, form the units for petrogenetic studies. 
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They give a correct picture of the fundamental chemical characteristic of 
the mineral or the rock, as a whole. They are valuable in interpreting the 
mode of origin of minerals, in making comparative studies of minerals and 
in the classification of rocks. They have the sanction of studies based on 
atomic chemistry. 


CLASSIFICATION OF IGNEOUS ROCKS 
Niggli Values.—For purposes of classification, Niggli has taken primarily 
into consideration : 
(1) alumina, 
(2) oxides of iron and magnesium, 
(3) lime, and 
(4) alkalies. 
These four items are represented by symbols which together with the 
mode of calculation are given below: 
(1) Determine the molecular numbers, 


Molecular numbers = weight percentage 
molecular weight 


(2) Add the molecular numbers of: 





1,000. 


Molecular  Niggli 
numbers symbols 


(a) Al,O3. Cr,O, and rare earths ae m’ al 
(b) Fe,O; X 2, FeO, MgO and MnO ... fm 
(c) CaO, SrO and BaO - = p Cc 
(d) Na,O, K,O and Li,O 
The total molecular numbers of 
al, fm, c and alk = (m' + n' + p’+ q') =x’ (say) 


‘he percentages are then determined thus: 


al = ™ x 100 
x 


’ 


fm = —, x 100 


c= +, X 100 


x 100 
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In addition to these four vaiues, other Niggli values are derived thus: 


si = = x 100 


Stated on per cent. { ti= mol. no. of TiO, x 100 


basis * 





_ mol. no. of P.O; 
x’ 


> x 100 


Stated as a ratio 


_ mol. no. of MgO 
‘ng fm 





_ mol. no. of KO 
| 


mg 
Thus the molecular values—al, fm, c, alk, si, ti, and p are expressed 
on a percentage basis, whereas k and mg are expressed as a ratio. 
qz, a derived value, known as quartz index is also used. 
qz = si — si’, where si’ = 100 + 4 alk, when alk < al. 
qz = si — (100 + 3 al + 1 alk) when alk > al. 


qz may be positive or negative and gives valuable clue of the minerals 
expected to be formed. 


On the basis of these values the eruptive rocks are classified into calc- 
alkali, soda-rich and potash-rich divisions. The divisions are subdivided 
into groups and magma types as mentioned below: 


Calc-alkali division, into 13 groups with a total of 64 types. 
Soda-rich division, into 15 groups with a total of 65 types. 
Potash-rich division, into 12 groups with a total of 45 types. 


One example of the calc-alkali division with the types is given below: 
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TABLE | 








Division: Group: 
Calc-alkali Hornblendite magma 





Types i fm c mm & mg 





al-hornblendite .. = 61 15 ap 0:7 





hornblendite ii io 60 20 aa 0:6 





hornblende peridotite -. & 74 14 its 0-7 





diallagite “i — 7 68 23 - 0:7 





websterite _ = a 1 80 18 ve 0:8 





Generally alk = 8 
al S 18-5 
fm > 50 
c— 9to 26 
si range ee al about 20 c<23 
si ‘ni al=10 c—15 to 26 fm<70 
si ‘a al about 12 c— 9 to 20 fm=70/69 
to 7-5 
si - al<10 c=20 fm nearly 70, alk 
practically nil 
si al<7°5 c— 9 to 22 fm >70, alk prac- 
tically nil 


TABLE IT 


Molecular numbers and Atom numbers 
Hornblende Norite—St. Thomas Mount, Madras 





Constituents e. Mol. No. x 1,000 Atom No. 





SiO, ae = ° 834 834 
ALO, .. . 114°5 229 
FeO, .. oe . 16°4 

FeO oe 7 , 218-5 

MnO 2% ve ah “a 

MgO me we ° 139-5 

CaO a nd , 141 

Na,O... we ° 50 

K,O a ae . 10 

TiO, ie e . 24 

P.O; we ee z 1 “0 

H,0O+ 


Total - na Analyst—Washington 















ilk 
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TABLE III 
Calculation of Niggli Values 











Niggli Molecular Niggli 
symbols numbers values 
al Al,O, + Cr,O, + Rare Earths 115 : 
wit © + 8 Se x 100 16-27 
fm FeO + Fe,0O, + MgO + MnO 391 , 
218-5+ 33 +139-5+4+ 0 “a tO. 
ec CaO +BaO +Sr0 141 
mam te +e a 141 707 x 100 19-95 
alk K,O +Na,O + Li,O _60 ; 
10 +a 4 © * - 60 707 x 100 8-48 
707 100-00 
, 834 
' 24 
ti ins ue Zh és 24 707 x 100 3-404 
K,O 10 
k K,0+Na,O ee ee ee ee 60 0 166 
MgO 139-5 
If the rock is an eruptive one it belongs to 
Group .. Hornblendite Magma } Refer values 
Type .. *Hornblendite ’ already given. 


If it is metamorphic it must have been derived from a calc-alumino- 
silicate rock with appreciable quantities of alkalies, say an alkali calc shale. 


CHEMICAL EQUILIBRIUM IN MAGMATIC MELTS 


Matter is said to be in equilibrium with its surroundings if it does not 
show any tendency to change. This condition of equilibrium is defined as 
a “system”. A system exists as such so long as the physical and chemical 
conditions remain unaltered. Such equilibrium conditions may. be obtained 
in homogeneous and heterogeneous bodies, In a homogeneous system 
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every part of the system has absolutely the same physical and chemical pro- 
perties as every other part. An unsaturated solution of sodium chloride 
or a silicate melt (magma) is an example of a homogeneous system. 


Heterogeneous Fquilibrium.—In over-saturated conditions sodium chlo- 
ride settles down from the solution. There comes a bounding surface between 
the solid and the solution. The physical and chemical properties of the 
solid are different from those of the solution. The two together constitute 
a heterogeneous system. This system is in equilibrium so long as the physical 
and chemical conditions remain the same. In the same way the moment 
the first mineral crystallizes from a silicate melt a new solid phase sets in, 
by which is meant a physically homogeneous but mechanically separable 
portion of a system. Thus the homogeneous system is changed to a hetero- 
geneous system. 


The mineral that crystallizes reacts with the liquid in which it is immersed, 
continuously or discontinuously, giving rise to additional minerals making 
the system more and more complex. 


A systematic study of the temperatures at which synthetic minerals are 
in equilibrium with silicate melts of known composition is one of the set 
programmes of the Geophysical Laboratory, Carnegie Institute cf Washington. 
The petrologic bearing of these studies is made out by equilibrium diagrams 


according to phase rule methods. The experimental facts have greatly added 
to the knowledge of the physico-chemical conditions relating to crystal 
differentiation in binary, ternary and quarternary systems. Bowen’s crystalli- 
zation differentiation theory has been generally accepted by distinguished 
petrologists all the world over. But the fact remains that the “crucible” 
experiments indicate only the conditions under thermodynamic equilibrium 
and fail to take into account “the kinetics of the processes involved in 
fractional crystallization ”’. 


Niggli is an advocate of Bowen’s differentiation theory. But he points 
out “ while all knowledge in this respect is unquestionably of the utmost 
value, the fact may not be overlooked that there must exist laws of a higher 
order, which govern the entire complex of phenomena ”’. 


Homogeneous Equilibrium.—Niggli searches for the cause, in the reactions 
that take place among the molecules while in the fluid state. An insight 
into the !aws that govern reactions in the melt state, i.e., under homogeneous 
equilibria conditions, would enable one to forecast what “* new phase must 
appear during crystallization, taking place under circumstances of any kind” 
—(Niggli). The birth of a mineral is determined by laws that operate in 
the pre-birth condition, 
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The laboratory experiments themselves are suggestive in this regard. 
The metasilicate clinoenstatite (MgSiO,) melts incongruently into the ortho- 
silicate Mg,SiO, at 1,557°C. There is a theoretical possibility of its further 
dissociation into MgO and SiO,. Does not this suggest that the ortho- 
silicate has very little tendency to dissociate, and lead to the hypothesis of 
its independent existence in the “melt”? The polysilicate orthoclase 
(KAISi,0,) melts incongruently into leucite (KAISi,O,) at 1,170°C. A still 
lower silicated mineral (KAISiO,) has been observed to exist in the study 
of the system K,O-AI,O,-SiO,. Again the dissociation of mullite (Al,SI,O,;) 
into corundum (Al,O,) and liquid at 1,810°C. leads to the possibility of 
the existence of Al,O, molecule in the melt. Hence can we not postulate the 
potential existence of K,SiO, and Na,SiO, in a melt, which on combining 
with Al,O, give rise respectively to KAISiO, and NaAISiO, ? 


The emphasis is on the fact that Mg,SiO,, KAISIO, and other ortho- 
silicates such as NaAlSiO, and Ca,SiO, are unstable while SiO, molecules 
are present in the melt. It is even possible that molecules such as K,SiO, 
and Na,SiO, give rise to KAISiO, and NaAlSiO, by reacting with Al,O,. 


Basis Molecules.—The fundamental compounds that exist in the magma 
are called basis molecules. There is no reason to suppose that these could 
be present in larger amount in a further state of dissociation. These bases 
on further reaction are considered to have given rise to new minerals. At 
this stage the following questions arise: 


(i) What are the characteristics of the basis molecules ? 
(ii) Are the reaction equations to be balanced by weight, or by 
volume ? 


Chemical reactions are based on the principle that each element has a 
definite combining weight or equivalent weight, and that the combinations 
occur in simple multiples of these quantities. The atomic weight is either 
identical with the combining weight or an integral multiple of it. 


Atomic weight = Equivalent weight x Valency. 


Hence the atomic weight is sometimes called the combining weight. 


Volume Unit.—But in petrological studies the first feature that strikes 
one on examining a rock megascopically or microscopically is the fabric 
of the rock, restricting the term “fabric” to define the sum-total of the 
textural and structural features of the mineral constituents. Petrogenesis 
seeks to explain the factors that have led to the development and propcrtion 
of minerals, as we see them in the rock. For most petrological work, the 
planemetric analysis is resorted to in order to give the weight percentage of 
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the oxides. In this analysis the volume percentage is obtained directly, 
Such being the case, it would be only appropriate to base reaction equations 
on a volume basis. 

It is now an established fact that the electro-positive character of the 
elements is the root cause of the chemical reactions leading to mineral forma- 
tion. The more electro-positive a metal is, the greater its affinity for a parti- 
cular non-metal or a radical. The atomic number denotes the number of 
protons or of positive electricity in the nucleus of the atom. If two minerals 
A and B combine to form C, it may be considered as due to the atomic num- 
ber of A having combined with that of B to give rise to that of C. C has 
an atomic number equal to the sum of the atomic numbers of A and B. 


Atom Numbers as Units.—\nstead of summing up the atomic numbers 
of the elements. there would be a justification in summing up the actual! 
number of electro-positive atoms themselves, if the atomic numbers of the 
elements are the same. But they are not of course the same. However, 
for calculation purposes, we could consider them to be the same, and treat 
them all alike if the atomic numbers of the reacting elements do not vary 
much. 


Thus the atomic number of magnesium is 12: that of silicon is 14. 


MgO +SiO,_ .. a .. MgSiO, 
= +O 7 - - 


14 26 


12 *e *e *. 12 
If the difference } is neglected. then 


MgO +SiO, .. oe .. MgSiO, 
: «A a nm oe 


It is interesting to note that the atomic numbers of the elements that 
enter into the composition of the rock-forming minerals show a range of 11 to 
26 only as seen below: 


Element Atomic No. 


Si a ‘ni —. 
Al - = ees: 
Fe La oe ea a 
Mg ‘ as ~~ 
Mn ae sa i Sa 
Ca s - a 
Na Eis ry «J a 
K “3 e a. ae 


l + 
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This range may be ignored as the actual variation in atomic numbers 
in all the elements taken together is 1 to 92. 


Mineral Formula.—As the sum of the electro-positive elements in a 
molecule wi!l give a number which denotes the maximum capacity of the 
molecule for reaction, it should be possible to state the formula of a mineral 
in a form, which besides indicating its chemical constitution would also 
connote its capacity for reaction. 


Let us take an example. Forsterite has the composition 2 MgOSiOg,. 
Its chemical constituents are Mg, Si and O. The number of electropositive 
atoms are 3(2 of Mg and 1 of Si). 


(By definition, a basis molecule has only 1 electro-positive charge. 
Hence the basis molecule forsterite is + Mg,SiO,; it is given the symbol Fo.) 


Thus the following table gives examples of the formule and symbols 
of a few basis molecules. 


TABLE [V 


Some Mineral Molecules 





No. of electro- Basis 
Mineral molecule positive atoms (Equivalent formula Symbol 
(Cations) unit) 





Mg,SiO, .. . 4 Mg,SiO, 
2 


KAISIO, .. 7 } KAISIO, 


NaAlSiO, .. ‘a i 4 NaAlSiO, 


CaAlAlO, .. i 4 CaAIAIO, 


SiO, a 2, ae SiO, 
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The formation of new minerals results thus: 


(a) By reaction of a basis with another basis. 
Example : 


Basis Basis 
Kaliophyllite Quartz 


Symbol... 7" K 


p Q 
Formula .. KAISiO, + 2 SiO, 
Niggli formula pee 3 Kp + 2Q 

(b) By reaction of a basis with a mineral. 

Example : 


Basis Mineral 
Kaliophyllite Orthoclase 


Symbol .. ‘< K Or 


p 
Formula KAISiO, + KAISi,O, 
Niggli formula ‘i 3 Kp + 5 Or 


(c) By reaction of a mineral with another mineral. 
Example : 
Mineral 


Mineral 
Wollastonite 


Enstatite 


Symbol .. a Wo En 


Formula Se CaSiO, + MgsiO, 
Niggli formula ‘ca 2 Wo + 2 En 


Mineral 
Orthoclase 


Or 
KAISi,O, 
5 Or 


Mineral 
Leucite 


Ec 
K,Al,Si,O,, 
8 Le 


Mineral 
Diopside 


The coefficients are obtained by summing up the electro-positive ele- 
ments. Thus in the example (5): 


3Kp=1K+1Al+1Si 
5Or=iK-+1Al+3Si 
8Lo=2K+2Al+ 4Si 


It may also be pointed out that since 
3Kp+ 2Q= 50Or 


30Kp+20Q= 500r 
60 Kp + 40Q = 100 Or 


and similarly in other cases. For there is a proportionate increase in the 


number of electro-positive atoms. 


In this mode of calculation, H, C, N, O and F—elements of the first 
horizontal series of the periodic table—are not counted. In sulphides and 
halides S and Cl are not counted, as they take the place of O. But these ar¢ 


counted as electro-positive in sulphates and chlorates, 
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This principle has a great advantage in metamorphic processes. For, 
§, HO and CO, need not be considered, as the addition of any in a forward 
reaction is balanced by the escape of the same in the reverse reaction. Hence 
such constituents are put within brackets. 


Examples.—({1) Muscovite and quartz, by metamorphism can change 
over to orthoclase, andalusite, and water. The reaction equations are: 


Muscovite + Quartz ~ Orthoclase + Andalusite + Water 
KAI,Si,0,,(0H), + SiO, —~ KAISi;0, + Al,SiO; + (H,O) 


7 Ms + 1Q > 5 Or + 3 Sil + (W) 


(2) Cordierite, orthoclase and water can change over to biotite, muscovite 
and quartz. The reaction equations are: 


Cordierite Orthoclase Water 
3 (Mg, Fe),Ai,Si;0,, + 8 KAISi,0, + 8 H,O 
33 Cord + 400r + (8 W) 


Biotite Muscovite Quartz 
—~—» 2K (Mg, Fe);AlSi,0,;,(OH), + 6 KAI,Si;0,,(OH), + 15 SiO, 


16 Bi + 42 Ms + 15Q 
Summary.—To summarise, chemical reactions can be conceived to take 
place in the melt condition. The melt is postulated to contain low silicated 
molecules such as Mg,SiO,, KAISiO,, etc. These plus basis molecules, by 
reaction, give rise to simple mineral molecules. The latter may combine 
among themselves, or with the basis molecules and lead to the development 


of other minerals. Even the growth of complex minerals can be explained 
by such reactions. 


The reaction equations are represented in a simple manner. As the 
electro-positive character of the elements is directly responsible for chemical 
reaction, the sum of the atomic numbers of the electro-positive atoms of the 
reactants and that of the resulting minerals could be balanced. But as there 
is scarcely any difference in the atomic numbers of the elements that go to 
build the rock-forming minerals, no great error comes in the equation if they 
are balanced by taking into account the sum of the electro-positive atoms 
on both sides of the equation. 


The basis and the mineral molecules are given symbols. Thus Q stands 
for quartz. Or means orthoclase. Reaction equations are made with co- 


efficients that denote the number of the electro-positive atoms in the consti- 
tuents, 
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Example : 
10 Or + 12 En + (2 W) > 16Bi+ 6Q. 


According to this line of thought the oxides which are more fundamental 
than even the basis molecules, from the reaction point of view, should have 
equivalent volumes. 


The following table shows that the oxides satisfy nearly this condition, 


TABLE V 
Rock-Forming Oxides 





Equivalent weight Equivalent volume 





SiO, _ .. 
1/2 Al,O, 
1/2 Fe,O, 
FeO on 
MnO 
MgO 
Cae ... 
1/2 Na,O 
1/2 K,O 
TiO, 


22-7 
12-9 
15-1 
11-8 
13-2 
11-3 
16-7 
12-8 
18-2 
18-9 


SAN UD 


~] wn 
BSSZRSESzSHS 
OHO Wonmnnwo-— 





The following table gives a list of ‘ Basis Molecules’ with their res- 
pective symbols. 
TABLE VI 


Basis Molecule Groups 





Basis Basis 
(Equivalent Symbol Equivalent (Equivalent Symbol Equivalent 
formula weight formula weight 
unit) unit) 





- Al Si O, 
Na Al Si 
'Ca Al Al 
|Mg Al Al 
| Fe Al Al 
1/3 \K K Si 
| Na Na Si 
|Ca Ca Si 
Mg Mg Si 
| Fe Fe Si 
\Fe Cr Cr 
1 Na 
Al Al 
1/2 be Fe 
Za Si 


es sIN 


1/3 Fe,Si 


PHRPYUSL“YSBSB 


aAA~ 
DAooprrnbhofrrfh=~ 





— 
Se Bee 
ore a 
wor 
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It is best to remind ourselves that the aim of the study is to trace the 
jine of development or growth of the rock as we have it for study. 


This involves: 
(1) The determination of the constituents either 
(a) by the planemetric analysis and the calculation to percentage 
weight, knowing the specific gravity; or 
(b) by actual chemical analysis. 
(2) The compounding of the oxides to determine the basis. 


(3) Compounding the basis to form the minerals we see which will 
be called katanorms. 


(4) Recognising the possibility of different lines of comoination which 
would yield different katanorms. 


{5) Finally synthesising them in such a way as to lead to the goal 
of modal percentage and drawing out inferences from the changes in the 
mineral constitution. That is, a suite of rocks may have a katanorm combi- 
nation A and may also have another katanorm B or “ katavariant”. If 
two rocks in this suite have their modal percentage corresponding to A and B, 
a consanguinitic relationship of one to the other is established. 


Or again, in a metamorphic suite, A, B and C may represent the kata- 
mezo and epizonal conditions. 


Petrogenetic studies involve the searching for the cause of the kata- 
normative variant corresponding to the modal mineral variant. 


It should be obvious that Niggli’s molecular norm—based on funda- 
mental values, is adapted to bring out the relationship between chemical 
composition and mineral! constituents of igneous and metamorphic rocks. 
The CIPW norm that is based on weight percentages of hypothetical norma- 
tive minerals sometimes leads to the determination of norms not at all 
present in the mode. Tom F. W. Barth observes: 


** The calculation of the so-called ‘ norm’ (CIPW norm) is a tool which, 
particularly for igneous rocks, has proved its great value. But in metamor- 
phic rocks, in which the interrelation of mineralogy and chemistry is of the 
utmost importance, the norm classification has failed. The reason is obvious: 
The normalization of the different metamorphic mineral facies encounters 
great difficulties in a norm based on weight percentages. The chemical 
relations of a rock in terms of weight percentages obscure the comprehensive 
view; the computations, moreover, become unnecessarily cumbersome. 
I propose, therefore, that calculation of the classical weight norm be altogether 
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discontinued. I regard it as obsolete and superseded partly by the Niggli 
‘molecular norm’, partly by the norm of the standard cell, as will be 
introduced in the present paper.” 


BASIS—RULES FOR THE CALCULATION 


From Atom Numbers 
1. Determine the molecuiar numbers of the constituent oxides 


Mint 18n. <- Weight percentage of the oxide 


~~ Molecular weight — x 1,000. 


2. Determine the number of atoms in each constituent from the mole- 


cular numbers as follows: 


The number of atoms are obtained by multiplying the molecular num- 
bers by 2 in the case of Al,O;. Fe,0,, Na,O, K,O and POs. 


3. The number of atoms which may hereafter be called ‘** atom numbers ” 
are equal to the molecular numbers in the case of SiO,, TiO,, FeO, MnO, 
NiO, MgO, and CaO. This also holds good for H,O. 


4. The atom numbers of: 
(a) MnO and NiO are added to FeO, 
(b) SrO and BaO to CaO, 
(c) Li,O to K,O, 
(d) Cr,O, to Al,O,. 


5. Appreciable quantities of P,O;, Cl, F and S are directly converted 
into basis molecules. 


P.O; will form the basis apatite .. (Cp) 
Cl will form the basis halite .. (AN 
F will form the basis fluorite .. (Fr) 
S_ will form the basis pyrite -+ ay 
6. Very small quantities of the constituents can be ignored. 
7. The number of atoms in each constituent is always referred to with 
a symbol, the first letter being in capitals thus: 


(Symbols of Atom Numbers) 


Constituents .. SiO, Al,O, Fe,0; FeOQ MgO CaO Na,O K,O TiO, 
Symbol ae Al Fe” Fe’ Mg Ca Ne K Ti 





Mg 
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Procedure 
A 


(i) Start with K and allot equivalent quantity of Al and S 
to make Kp which has the formula 1/3 KAI SiO,. Thus, 


xK + xAl + xSi = 3xKp = 3 xf APO 


(ii) If K > Al, the balance of K after (i) is converted into the 
basis 1/3 K,SiO, = Ks, thus, 


Yan 3 os 1] 
YK + 5 Si = 5 Ks = 5 3 
(iii) If Al remains after the saturation of K, then proceed with 
the calculation of basis Ne with Na. 


B 


(i) xNa + xAl+ xSi = 3xNe = 3x [MA SO*1. 


(ii) If Na is in excess of Al, the balance of Na is converted 
into the basis Ns. 


: 3 a Na, SiO. 
yNa + 5 Si= 5 Ns = a ‘| 


(iii) On the other hand, if Al is in excess of Na, then turn to Ca 
Cc 


(i) ) Ca + xh 5 xCal = 5 [| z 7 


(ii) If Ca remains over after (i) above, then form Cs. 


Y om 3 5 = 3 y [22804], 
yCa + 9) si = 7 VCs -_ y =| 


(iii) If on the other hand, Al remains in excess, then turn to Mg 


D 
(i) The basis Sp is formed thus 


y 3 = boon ; 
5 Mg + yAl= 59Sp = 5 3 ip 
(ii) If after Sp is formed, there is yet a balance of Al, then 
the saturation of Fe is taken to form the basis hercynite (Hz). 


(iii) If on the other hand there is a balance of Mg proceed 
as directed in (G. i), after going through E and F. 
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E 


2 — 3 


(ii) If still Al is in excess, corundum is formed with the 
balance say, w 


w. Al= w.C = 


(i) 5 te’ + a ee [= 3). 


w. Al,O; 
| a 


(iii) If Fe” is in excess, proceed as directed under (G. ii) after 
taking the steps F and (G. i). 





F 
If Fe’”’ is present, then convert it to form the silicate Fs thus 


xFe” + 3 Si= 5 xFs=3} x [Fee]. 


G 


| (i) The balance of Mg in (D.iii) is converted into the basis 
forsterite (Fo) 


x.Mg + 5 Si =5 x. Rom x[ Mee?.1. 


2 


(ii) After the above step, the balance of Fe” in E is converted 
into Fa as follows: 


«Ro 2 _ 3 [Fe,"SiO,1. 
xFe + 3 Si= 5 xFa = 5 x [2] 


Of the minor constituents, 
(i) TiO, is considered as basis rutile (Ru). 


(ii) ZrO, is combined with an equivalent amount of SiO, to 
form the basis 1/2 ZrSiO,, Zircon (Z). 


xZrO, + xSiO, = Ysa A a 2x{ A504) L 


~ 


(iii) Lastly, the balance of SiO, is considered as basis Q. The 
other minor constituents are disposed of thus: 


xP,O, to form apatite (Cp) 


xC! to form halite (H1) 
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xF to form fluorite (Fr) an SF 


xS to form pyrite (Pr) = xPr 


xCr to form chromite (Cr) = ; xCr 


It would therefore be possible to arrive at the following groups 
of combinations. 


TABLE VII 
Basis Groups 


(in relation ‘to atom numbers) 





Alkali excess 


K-+Na K+Na<Al 
+2 Ca also Alc K+Nat2 Ca 4). K+Na+2 Ca 


K<Al Al. K+Na+2Ca>Al +2 Fe’+2Mg “~ +2 Fe’+2 Mg 


Al also Al>K-+Na-+2 Ca 
K+Na>Al 





2 5 





Kp Kp 
Ne Ne 

Ns Cal 

Cs Sp 

Fs Fo or Hz 
Fo Fa 

Fa Fs 

Ru Ru 


Q 


Sum up the atom numbers of all constituents. 





Determine the percentage of each basis molecule as in 


: x 100, etc., where 


Example.—Hornblende Norite—St. Thomas Mount, Madras. (Analyst 
—Washington). 


Calculation of the basis by distribution of atom numbers. 
A glance at the atom numbers shows that 

K + Na < Al 

K + Na + 2Ca > Al 


Hence, the basis of Group 4 will be formed, 
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I From Niggli Values 


The basis molecules have been derived from the atom numbers. It is 
possible to calculate them directly from Niggli values if they are already 
available. 


This involves setting forth the relationship between molecular numbers 
used in Niggli values and the atom numbers that form the units for the 
calculation of the basis. 


The following table gives the relationship: 
TABLE IX 
The Atom—The Molecule—Equations 





II Ill IV 





Symbols 
Constituents Relationship Reference 


Appendix B 
Atoms Molecules 





Si Si =si 

Al Al 2 al 

Fe’’’ 

Fe’ \ fm Fe =fm.(l—mg) 
Mg mg Mg = fm (mg) 
Ca c Ca 4 

Na as Na 


K k K 
Ti Ti 


alk (l—k) 
alk-k 


oo MNO 


oo 


Hou te dl 





2 = si+2 al+fm-+c-+2 alk+ti 


But al-+-fm-+c--alk = ™ oat PT y 100 (see p. 2) 
= 100 





+ = si+al+alk+ti+100 
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TABLE X 
Basis—Nigegli Values, Equations 











Basis Niggli values Reference 

Appendix B "7 

KAISiIO, .. .. Kp 6 alk-k (5) 

NaAISiO, .. .. Ne 6 alk (I—k) (6) 
CaAlSiO, .. — 3 (al—alk) (7) B 

CaCaSiO, .. . ¢s ; (c—al-talk) (8) 
1/3 i 
| ) 3 al 
MgMgSi0, .. .. Fo 5 fm-mg (9) fn 
¢ 
al 

| FeFesiO, Fa 

SiO, a si—(50-+2 alk—al) (11) " 
TiO, ~ << ti ti 





= = 100+si+alk+al-+ti 














Niggli’s Principles of Igneous Petrogenesis 


Calculation of the Basis from the Niggli Values 





Basis percentage 





Niggli values Basis from Niggli value equations =(2) x 100+246 
=(2)x -406 
(1) =(3) 
Kp =6alk-k 
= 6x 8-48 x +166 = 8-450 8-45x-406= 3-431 
¢ =118-00 Ne =6alk(I—k) 
al = 16-27 = 6x 8-48 x -834 = 42-44 42-44 -406 = 17-230 
fm = 55:30 Cal = 3x (al—alk) 
“oe *~ = 3x (16-27—-8-48) = 23-37 23-37x -406 = 9-488 
k= 8-48 
oe G& « ; (c—al-talk) = 18-24  18-24x -406 = 7-404 
mg= 0°30 Fo = } fm-mg 
i= 3-4 ™ ; x 55-30 0-36 = 29:86  29-86x -406 = 12-130 
3 
Fs+Fa = 5 fm (l1—mg) 
=3 x 55-30% +64 = 53-09 53-09 -406 = 21-550 
Ru =ti = 3:40  3-40x-406= 1-380 


= si—(50+2 alk—al) 


= 118—(50+16-96—16-27) = 67-31 








246-160 





67-31 x -406 = 27-330 


99 -943 
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KATANORMATIVE MINERALS 


The minerals of the basis groups are to be considered as fundamental] 
units which have to combine with silica to form the saturated minerals. The 
minerals formed from the basis in the first instance, according to the prescribed 
rules mentioned below are called katanormative minerals or katanorms. 


The katanorms thus formed may or may not be present in the rock. 
If they do not give us a picture of the modal percentage of the minerals, the 
katanorms have to be split or synthesised in such a manner that the modal 
percentage of minerals results out of the reaction. In this process there will 
be stages of growth, each stage representing a particular katanorm variant. 
The aim, however, is to get at the variant corresponding to the mode. 


Katanorm— Rules for the Calculation 








Bases Reaction Katanorm 





3 Ns 1. Ns is combined with an equal quan- Fns 
(Na, Na, SiO,) tity of Fs to form: 2Fns. (Fns) SiO,NaFe’’ 


3 Ks Ks. Do. 2 Fks. Fks 
(K, K, SiO,) (Fks) SiO,KFe”’ 


(a) 3 Fs = Fe,SiO, 
= Fe,0,+Si0, 
=2Hm+1Q Hm—Hemiatite 
(Fe,0,) 
(b) Mt and Q will result when two 
units of Fs interact with one 
unit of Fa as follows: 
3 Fs 6 Fs+3 Fa = Fe,Si,0O,)+Fe,SiO, 
(Fe,SiO,) = Fe,Si,0,,4 
=> z Fe,0,+3 SiO, f 
=6 Mt+3Q Mt—Magnetite 
Therefore 2 Fs+Fa=2 Mt+1Q Fe,O, 


(c) If more than two units of Fs 
are available for each unit of 
Fa, then, both magnetite (Mt) 
and hematite (Hm) will result 
as follows: 

Thus if x units of Fs are 
left over from 1, let 2 y units 
of these combine with y units 
of Fa to form 2 y units of Mt 
and y units of Q. The balance 
of (x—2 y) units of Fs will 
form 4(x—2y) of Hm and 
4 (x—2 y) of Q. We therefore 
have 








—_—_—————— 
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Reaction 


Katanorm 





3 Cal 
(CaAl,0,) 


3 Cs 
(Ca,SiO,) 


Fks 


Fns 


3 Ns 
(Na,SiO,) 


3 Ks 


(K,SiO,) 


3Kp 
(KAISiO,) 


x-Fs+y:Fa 
= (x—2 y) Fs+y (2 Fs+Fa) 
= 2 y-Mt+4 (x—2 y) Hm 
+Ly+4 (x—2 y)] Q 
where x>2 y [vide (a) & (b) above] 


The liberated Q must be 
added to the originally present 
Q and distributed. 


If Cal is present, anorthite is ob- 
tained. 


xCal+% xQ= 3 x-An 
If Q is present after (3), wollastonite 
is formed. 
x-Cs-+4x-Q= 3 x-Wo 
If Q is present after (4), potassic 


acmite (K-Ac) and Sodic acmite 
(Ac) are formed. 


x-Fks+ 3 Q= f x-K-Ac 


a, Se 
x:Fns+ 3 Q= 3% Ac 
If during the formation of Fks and 


Fns, Ks, or Ns are left over, then, 
x-Ns +42xQ= 4 x-Ns’ 


x:Ks +4xQ= ; x*Ks’ 


If after all these operations Q still 
remains, then orthoclase is form- 
ed. 

x-Kp+3 xQ= 3 x-Or or 


Kp may have to be given to a 
balance of Q, say y, in which case 


3 5 
YQ+ 5yKp= 5 y-Or 


Mt & Hm 
Magnetite & Hematite 


An 
Anorthite 


Wo 
Wollastonite 


Ns’ 
Sodium-di silicate 
(Na,Si,O;) 


Similarly Ks’ 


Orthoclase 
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Reaction Katanorm 





8. For a further excess of Q albite is 
formed. 


3 Ne 3 5 Ab 
(NaAlSiO,) yYQt 5 y'Ne= 5 Ab Albite 


9. If Sp and Hz bases are present and 
if there is still Q left over, cordie- 
rite (cord) is formed. 


3 Sp 6 11 Cord 
(MgAl,O,) yQ+ 5¥°Sp= 5 y-Cord Cordierite 


3 Hz il 
(FeAl,O,) i Q+8 y'Hz= = y:Fe—Cord 


10. The basis molecules that are left 
over are Fo, Fa and C. To the 
3 Fo extent possible, En 
(Mg,SiO,) (i) Fo is converted into enstatite Enstatite 
(En). Thus 
y'Q+3 y-Fo=4 y-En 
3 Fa (ii) Fa is converted into hyper- Hy 
(Fe,SiO,) sthene Hypersthene 
y'Q+3 y:-Fa=4 y-Hy and 
2c (iii) C is converted into sillima- 
(Al,O,) nite (Sil) Sil 
y'Q42 y-C=3 y-Sil Sillimanite 





Inasmuch as the basis is calculated to a total of 100, the new assembly 
of minerals will also give a total of 100. 


As has been already said, the ultimate object is to get at such minerals 
as are actually present in the rock. In the transitional stages, the kata- 
normative minerals may have to pass through hypothetical katavariants. 


In special cases, the katamolecular norms are calculated with reference 
to specia! theoretical mineral compounds. 
p 


Tke mineral compounds may be made up of a basis molecule and a 
mineral or two minerals. The following examples are illustrative. 


3Kp +5Or = 8Le (Leucite) 

3 Cal + 2Wo = 5 Ge (Gehlenite) 
3Cal+1Q =4Ts (Tschermak’s silicate) 
3Ne +8En +5Ab+6Fo 
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EXAMPLE 
Hornblende Norite 


Calculation of the Katanorm 








+ Ne Cal Cs Fo Fs Fa Ru Q 
.. 3-431 17-230 9-488 7-404 12-130 2-832 18-718 1-380 27-330 


Total = 99-95 











Calculation Katanorm UsedQ Released Q 





2Fs+1 Fa 


-2 Mt-+1Q 
7-832 Fs+1-416 Fa = 2 
5 
5 


32 Mt-+1-416Q Mt 


Or 


8 2-832 
718 Or Or 


3Kp+2Q 
3-431 Kp +2-287 Q 


: 5-718 2-287 


3Ne+2Q Ab 


Ab = 28-716 11-487 


3Cal+2Q 
9-488 Cal +6 -325 Q 


3Cs+1Q 
7-40 Cs-+2-46 Q 


5 

17-230 Ne+11-487 Q=28-716 Ab 
5 
5 


An 
-813 An 


4Wo 
9-86 Wo 


1 


An = 15-813 6-325 


Wo= 9-860 2-460 





22-559 





mbly Total Q available = 27-330 — 22-559 + 1-416 = 6-187. 
3Fo+1Q=4En 
3Fa+1Q=4Hy for, 
1Q=4(En + Hy) 


6-187 Q = 24-718 (En + Hy) 


1erals 
kata- 
nts. 


, Fo . 
En proportional to + Po (available) 


o 12-130 oe 
= 19-130 + 17-302 * 24°78 = 10-21 


En = 10-21 
Hy = 14-54 
Fo used for En = 10-21 x 3 = 7-56 
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Fa used for Hy = 14-54 x 3 
= 10-905 
Fo = 12-130—7-650 
= 4-48 
Fa = 18-718—(10-905 + 1-416 to: 
= 6:40 
1Ru+2Wo =3Tn 
1:38 Ru + 2-76 Wo = 4:-14Tn 
Tn =4-14. 
Balance of Wo = 9-86—2-76 = 7:10 
Wo =7:10. 


Theoretical Katanorm of Hornblende Norite (1) 





Katanorm Or Ab An Wo En Hy Fo Fa Tn Mt Total 
Percentage 5:72 28-72 15-81 7-10 10-21 14:54 4-48 6-40 4-14 2-83 99.95 





The katanorm thus derived according to prescribed rules does not 
reflect to any great extent the mode of the rock, which is composed of 
andesine, diopside, ortho-pyroxenes, hornblende, and iron ores. Further, 
sphene, wollastonite and olivine are not present in the rock. 


So an attempt must be made to evolve a katavariant which corres- 
ponds to that of the mode. [!n doing so, it is always desirable to take into 
account the actual analysis of the mineral and its volume percentage if 
available; a micrometric analysis indicates 20-25% hornblende. 


As there is no sphene, in this case, Tn can be split back. 
4-14 Tn = 1-38 Ru + 2-76 Wo 


and the katanorm, which is really a katavariant may be rewritten as 
follows :— 
Katavariant (2) 





Katanorm Or Ab An Wo_- En Hy Fo Fa Mt_ Ru Total 
Percentage 5-72 28-72 15-81 9-86 10-21 14-54 4-48 6-40 2-83 1-38 99-95 


eee 
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Diopside can be formed thus: 
2 Wo + 2 (En + Hy) = 4Di 
9-86 Wo + 9-86(En + Hy) = 19-72 Di 


, , 10-21 
En in Di = 94-75 x 9-86 


= 4-07 
Hy in Di = 5-79 
Deducting Wo, En and Hy used, and including Di obtained, tne 
katanorm will be as follows: 


Katavariant (3) 





Katanorm = Or Ab An Di En Hy Fo Fa Mt Ru _ Total 
Percentage 5-72 28-72 15°81 19-72 6-14 8-75 4:48 6-40 2-83 1-30 99-95 





There is no Q to form Ho which is 20 to 25% in the mode. 
Ho is obtained as follows: 


8 Di+ 6(En+ Hy) +1Q+(1W) =15Ho 
Q required for about 20% Ho = 1-35. 


This 1-35 Q can be liberated from Or, Ab and An, proportionately 
according to their availability. 


TABLE XII 
Ratio of Released Quartz from Feldspars 





Available from Quartz released 
(1) (2)=(1) x 1-35+50-25 





Or 5-72 0-154 
Ab 28-72 0-771 
An 15-81 0-425 





* Total 50-25 1-350 
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TABLE XIII 
Ratio of Released Kp, Ne and Cal 








Quartz Added Constituents 
released constituents to be deducted Katanorm  Katavariant 
from to Ho from (5) =(5)—(4) 








(1) (3) (4) (5) (6) 





Or = 0-154 x >= =0-231 Kp Kp =0-231 Or =0°385 Or = 5-72 Or= 5-3 


Ab = 0-772 x >= =1-156Ne Ne =1-156 Ab=1-:927 Ab=28-72 Ab=26-79 


An = 0-425 x 5 =0-638 Cal Cal=0-638 An=1-063 An=15-81 An= 14-75 


NIW NIW NI Ww 





Since, 

8 Di+ 6 (En + Hy) + 1Q+ (1 W) = 15 Ho, 

10-8 Di+ 8-1(En + Hy) + 1-35 Q + (1:35 W) = 20-25 Ho. 
8-1 (En + Hy) is made up of 3-34 En + 4-76 Hy. 


1 available En. | an : 
En = svailable (En + Hy) * °°! = ia-g9 * ©"! 


= 3-34 
Hy = 4-76. 


Katavariant (4) 








Katanorm Or Ab An Di En Hy Ho Fo Fa Mt Ru Total 
Percentage 5-33 26-79 14-75 8-92 2-80 3-99 22-28 4-48 6-40 2-83 1-38 999% 





At this stage, it is obvious that as there ‘s no olivine in the rock, Fo and 
Fa should be converted into the meta-silicates En and Hy. 


The feldspars are again indented for SiO, required. The Kp, Ne and 
Cal liberated will be fixed up again in hornblende (Ho). 


4-48 Fo+ 1-49Q = 5-97 En 
6-40 Fa + 2-13 Q = 8-53 Hy 
Total Q = 1-49+ 2-13 
= 3-62 
3-62 Q must be taken proportionately from the available Or, Ab and An, 
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TABLE XIV 


Ratio of further released quartz from feldspars 





Available from Quartz released 
(1) (2)=(1) x 3 -62+46 -87 











TABLE XV 
Ratio of Further Released Kp, Ne and Cal 





Added Constituents 
Basis constituents tobe deducted Katanorm  Katavariant 
to Ho from (5) 6=(5)—(4) 
=(1)-+(3) 





(1) (2) (3) (4) (5) (6) 





Or=0-41 x5>=0°61 Kp Kp =0-6l Or =1:02 Or= 5-33 Or= 4-31 


Ab=2:07 x =3-11Ne Ne=3-11 Ab=5-18 Ab=26-79 Ab=21-61 


Ni w NIW NIw 
I 


An=1:14 x 5>=1-71Cal Cal=1-71 An=2°85 An=14-75 An=11-90 





Therefore Ho (fina!) = 22-28 + 5-43 = 27-71. 


The final katavariant values are 


Katavariant (5) 





Constituents .. Or Ab An Di En Hy Ho Mt Ru _ Total 
Percentage .. 4°31 21-61 11°90 8:92 8°87 12°52 27-71 2-83 1-38 100-05 





It may here be remarked that the necessary Q for the formation of Ho 
may be released by forming Bi. But in such a process, calculation shows 
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that there is no adequate Q for the formation of about 25 per cent. horn- 


blende; further, it may not perhaps be quite appropriate to form Bi first 
and later Ho. 


CALCULATION OF THE BASIS OF HORNBLENDE 


The hornblende derived has a volume percentage of 27-11 and its 
composition is shown below: 





Constituents - Kp Ne Cal Di En Hy Qs‘ Total 
As a ratio .. 0°84 4-27 2-35 10-80 3-34 4-76 1-35 27°71 





The constituents Di, En and Hy, have to be split into the bases Cs, Fo, 
Fa and Q as follows: 


10-8 Di = 5-40 Wo + 2-48 En + 2-92 Hy. 
(The ratio En: Hy is taken from the primary katanorm.) 


Then, the composition wil] be 





Constituents .. Kp Ne Cal Wo- En Hy Q ‘Total 
As a ratio .. 0°84 4:27 2:35 5:4 3-34+ 4-76+ 1-35 27°71 
2-48 2-92 
5-82 7-68 





5:40 Wo = 4:05 Cs + 1:35Q 

5-82 En = 4:37 Fo+ 1:45Q 

7:68 Hy = 5:76 Fa+ 1-92Q 
Thus the bases are: 





Constituents .. .. Kp Ne Cal Cs Fo Fa Q_ Total 
Asaratio .. .. 0°84 4-27 2-35 4:05 4-37 5-76 6-07 27-71 





Inclusions of magnetite have not been taken into account in the bases. 
Mt = 2-83 


Allowing approximately 1-4 Mt for the hornblende as inclusion, 
1-40 Fs + 0-70 Fa = 1-40 Mt + 0-70Q 
Thus, we add in the basis of Ho 1-40 Fs and 0-70 Fa and deduct 0-70 Q. 
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TABLE XVI 
Hornblende—Basis 


Weight per cent. 





Constituents As a ratio sc x 3-435) C.R.’s analysis 





(3) (4) 





2-89 
14-67 
8-11 
13-91 
4-81 
14-97 
22-19 26: 
18-45 23-60 


— 
NWN OON 


essyak 


— 





97-90, 
Total .. 100-00 100-00 (approx.) 





It is very interesting that the calculated bases of column 3 are in close 
agreement with the bases given by Rajagopalan, vide in column 4, of horn- 
blende from the same area. 

CONCLUSION 


The principles of Niggli’s method outlined for the study of the origin 
and growth of igneous minerals and rocks is also applicable to the study of 
the metamorphic minerals and rocks. The genetic relationship between, say, 
a suite of calc gneisses with calcareous sediments or of aluminous hornfels 
with aluminous sediments can be cogently and intelligently established from 
the beginning to the end. A picture of the stages of metamorphic evolution 
could be drawn. The kata, meso and epizonal norms of a suite come to 
have great significance. 


His methods are coming to be used more and more by petrologists for 
their interpretative studies. There is life in the variation diagrams pre- 
pared by using Niggli’s fundamental units. He writes, ‘‘ The resistance often 
met with in petrographic circles against a scientific and critical revision of 
analytical data probably has various causes. Beside the habit of merely 
taking over the analyst’s figures, the fact that a choice must be made between 
several equally good methods of recalculating analysis-values may have 
produced a tendency to use none at all.” 
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In Appendix A are given examples of reaction equations of simple and 
complex minerals. Some of the more important paragenetic mineral 
assemblages, under epi, meso and katazonal environment are also given. 


To conclude, Niggli’s treatment of the subject has a unique simplicity 
and beauty about it, as the method is directed, to understand the probable 
but hidden cause of the origin of minerals and rocks, before they were 
ushered out of the melt. 


SUMMARY 


Modern conceptions cn the processes of crystallization of igneous rocks 
and the processes of recrystallization of metamorphic rocks, have resulted 
primarily from a knowledge of equilibrium of investigated silicate systems, 
which mean heterogeneous equilibria conditions. But Niggli considers 


that homogeneous equilibria conditions in the fluid state control mineral 
formation. 


A critical insight into the factors that govern the formation of minerals 
is possible by adopting atomic and molecular units as demanded by studies 


in atomic chemistry. Niggli’s classification of rocks is based on Molecular 
Values. 


Tt has become the accepted practice to use proporticns of the molecules 


such as SiO,, Al,Os, K,O, etc., in petrological calculations. But Niggli 
recommends the use of proportions of the cations (silicon and metals). The 
calculations are made to a total of 100. 


Reaction equations are made with coefficients that denote the number 
of the electro-positive elements. The molecular amount of a mineral thus 
becomes proportional to the sum of the cations. 


Niggli’s view is that the modal minerals are evolved by reactions in the 
melt state and that therefore petrogenetic problems should be directed to 
get at the mode by a series of reactions. The various stages of the reaction 
correspond to katavariants, which finally end in a katanorm that can be 
equated to the mode. Igneous petrogenesis needs drastic changes. The 
petrographic data must be critically and scientifically studied with a view to 
the discovery and the expression of natural laws, reliance being placed on 
atomic chemistry. 
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APPENDIX A 
IMPORTANT ROCK-FORMING MINERALS 
Polysilicates 


Orthoclase és 3Kp+ 2Q=50Or 
Albite * 3Ne+ 2Q=5 Ab 
Anorthite cr 3Cal+ 2Q=5An 


Metasilicates 


Leucite a 3Kp+1Q=4Le 
Wollastonite .. 3Cs+1Q=4Wo 
Hypersthene ia 3Fa+1Q=4Hy 
Enstatite - 3Fo+ 1Q=4En 
Diopside ~ i 2 Wo + 2En = 4Di 
Hedenbergite .. 2 Wo + 2 Hy = 4 Hed 


Acmite mS 
ponte .. Ac 3Ns + 2Hm+3Q=84Ac 


Jadeite oe ae 3Ne+1Q=43Jd 


Augite .. Aug 4Di+ 3Sp+ 1Q=8 Aug 
(special composition) 
Hornblende ~- Ho 8 Di+ 6En+ 1Q+ (1 W)=15Ho 
Cordierite .. Cord 6Sp+5Q= 11 Cord 
Fe-Cord 6 Hz + 5 Q = 11 Fe-Cord 


Orthosilicates 


Cancrinite .. Canc 18Ne+2Cc= 20Canc 
Sodalite .- Sod 18 Ne + 2 HI = 20 Sod 
Hauynite .. Hau 18 Ne + 4A = 22 Hau 

CaSO, anhydrite.. A 

Nosean .. Nos 18 Ne + 3 Th = 21 Nos 

Na,SO, thenardite Th 
Grossularite ax ae 5 An+ 3Cs = 8Gro 

Pyrope oo yp 11 Cord + 8En= 16Pyp+3Q 
Almandite .. Alm 11 Fe-Cord + 8 Hy = 16Alm+ 3Q 
Andradite .. Andr 6Wo+2Hm= 8 Andr 
Forsterite .. Fo Basis 
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Fayalite io Basis 
Larnite « Se Basis 
Monticellite .- Mont 3Cs+3Fo= 6 Mont 
Meionite .. Me 15 An+ 1 Cc = 16 Me 
Marialite io oan 15 Ab + 1 HI = 16 Ma 
3NaAlSi;0,—NaC! 
3NaAlSi,0, 

4Na,SO, re 30 Ab + 3 Th = 33 Ma 
3NaAlSi,O, 

4.Na,CO; il 30 Ab + 2 Nc = 32 Ma 
Gehlinite - & 3 Cal + 2 Wo = 5Ge 
Akermanite . AR 3Cs + 2En = 5 Ak 
Vesuvianite oo Vee 18 (Cs, Fo) + 6Sil+ 1Q+ (2 W) 


Ca, Al, (Mg, Fe), = 25 Ves 

Si,O,, (OH, F), 
anata 
Andalusite 
Kyanite 
Zoisite co aa 9Cal+ 2Wo+5Q+ (1 W) = 16Zo 
Epidote -- Ep(*= 

Fe-Ep 6Cal+4Wo+ 2Hm+ 4Q + (1W) 
= 16 Ep (J 


Sil 2C+1Q=3Sil 


Analcite -» Ame 4Jd+ (1 W) =4Anc 
Staurolite .. Staur 3Hz+ 3Sil+1Q+ (1 W) =7Staur 
Muscovite .. Ms 10 Or + 4C + (2 W) = 14 Ms 
Biotite .. Fe-Bi 6Kp+9Fa+1Q-+(2W)= I16Bi 
Phiogopite .. Mg-Bi= 

Pgp 6Kp + 9 Fo+ 1Q+ (2 W) = 16 Pgp 
Siderophyllite .. Sdro 6 Kp + 10 Hy + 3 Sil + (2 W) 
(H,K.Fe,;Al,Si;O.,) . = 16Sdro + 3Q 


Chloritoid 
seal Ottr 3Hz+1Q+(2W)=4Oter 


Serpentine -. Semp 3 Fo + 2 En + (2 W) = 5Serp 
Talc aa 6En+ 1Q+ (1 W)=7Tc 

~ Kaolin . 3 Sil+ 1Q+(2W)=4Ka 
Titanite oe Tae 1 Ru + 2 Wo = 3 Tit 
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Tschermak Silicates a 


CaAlSiO, .. Ts 3Cal + 1Q=4Ts (49) 
MgAl.SiO, .. Mg-Ts 3Sp+1Q=4 Mg-Ts (50) 


EXAMPLES OF COMPLEX REACTIONS 
Paragenetic Associations 


A. Epizonal Minerals 


1. Serpentinization 6 Fo + (2 W) + (CO,) = 5 Serp+1 Mgs (51) 
Magnesite .. Mgs 9Fo+ 1Q+ (4W) = 10Serp (52) 
Periclase cn ae 15 Fo + (4 W) = 10Serp+ 4Ps+1Q 


removed in solution (53) 


2. Steatization .. 10 Serp + (3 CO,) = 7 Tc + 3 Mgs 


+ (2 W) (54) 
(Grammatite) Gram 
Tremolite 15 Gram+(4 CO,)=4 Dol + 7 Tc +4 Q (55) 
Dolomite .. Dol 





3. Steatization— 
Carbonatization 10 Serp + (6 CO,) = 3 Mgs + 7 Tc 
+ (3 CO.) + (2 W) 
or 6 Mgs + 4Q + (4 W) (56) 


4. Kaolinization 


Kaolin .. Kaol 
Potas. Carbonate Kc 30 Or + (6 W) + (3 CO.) 
= 12 Kaol + 6Kce+ 12Q (57) 
5. Albitization .. 5An+2Nc+4Q=10Ab+1Cc (58) 


Sodium Carbonate Nec 


Epi-Mesozonal Minerals 


6. Decarbonization 6 Dol + 4Q+ (W) = 7Tc+ 3Cc 
(Progressive) +(3CO,) (59) 
14 Tc + 3 Cc = 15 Ho + 2 Dol 
+ (CO,) + (H,0) 


!) 
2) 
3) 


t) 
5) 


8) 


9) 


0) 








Carbonization .. 


Anthophyllite 


Tremolite 
Almandite 


Staurolite 
Muscovite 


Forsterite 
Cordierite 


Sapphirinization 


Sapphirine 
Sillimanite 


Enstatite 


Pyrope 
Spinel 


Sapphirinization 


Phlogopitization 
Phlogopite 
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3Cc + 8Bi+6Q=12Di+50r 
+ (1 W) + 3CO,) (61) 


Mesozonal Minerals 


Anth 


Ho 
Alm 


Staur 


Fo 
Cord 
Sphr 
Sil 
En 


Pyr 
Sp 


45 Anth + 56 Di+ 4Q + (4W) 
= 105 Ho (62) 


8 Hy + 11 Fe-Cord =16Alm+3Q (63) 
8 Fe-Bi + 6 Sil-+ 1Q=7Ms-+ 8Alm (64) 
8 Alm + 15 Sil + (1 W)=21 Staur+2Q (65) 
8 Fe-Bi + 21 Sil + (W) = 7 Ms 


+ 21 Staur+1Q (66) 
11 Fe-Cord + 2 Hm + (W) = 7 Staur 
+3Mt+3Q (67) 


Katazonal Minerals 


4Dol + 1Q=2Cc+ 3 Fo+(2C0,) (68) 
42 Ms + 16 Bi+ 15Q = 400r 
+ 33 Cord + (8 W) (69) 


3Fo+ 12En+ 16C+1Q 
= 11 cord + 21 Sphr (70) 
10 An+ 12Fo+4Di+3Q 
= 8 Gro+ 18(En+ Hy)+3Sil (71) 
10 An+ 18 Fo + 4Di+3Q 
= 8Gro+ 3Sp+ 24En (72) 
10 An + 16 En = 16 Pyp + 8 Di + 2Q (73) 
5An+ 4Di= 3Sp+4Wo+2Q (74) 


15 An + 28 En + 12 Wo + 36Sp 
+5Q-+ (2 W) = 30 Anth + 42 Sphr 
+ 24 Gro (75) 


(Mg-Bi) 18 Sp + 20En+ 4Q + 4KO, + (4W) 


Pyp 


= 14 Sphr + 32. Mg-Bi (76) 


T. N. MuTHUSWAMI 
APPENDIX B 


= FeO + 2 Fe,0O, + MgO 
Fe = FeO + 2Fe,0, 
Fe =fm—MgO 


MgO 
=~ 


Fe = fm—fm-mg 
Fe = fm(1—mg) 


_ MgO _Mg 
ea 


Mg = fm-mg 


_ po 
= Na,O + KO 


Na,O = alk—K,O 
K,O= k (K,O + Na,O) 
= k-alk 
Na,O = alk—k-alk 
= alk (1—k) 
Na = 2alk(l—k) 


K 


{Na,O — molecular number to be multiplied by 
2 to get atom numbers. alk (I—k) is 
molecular number.]} (3) 


K,0 =k-alk 
K = 2alk-k 


(K,O0— molecular numbers to be multiplied 
by 2 to get atom numbers. k-alk value is 
in molecular numbers.) (4) 


K = 2 alk-k 
Kp =3K= 6alk-k (5) 








(1) 


(2) 


(3) 


(4) 


(9) 
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Na = 2alk(1—k) (see 3) 
Ne =3Na 
Ne = 6alk(I—k) 


Atom numbers of Al,O, are distributed to CaO, 
only after satisfying K,O and Na,O. 


Available Al = 2 (al — alk) 


3 
= 3 (al — alk) 


Cs is formed from CaO, left over after CaO 
allotment to Cal 


Ca allotted to Cal = al—alk 
Ca balance for Cs = (c—al + alk) 


Cs = 5 (e—al + alk) 
3 
Fo = 5 Mg atoms 


= ; fm-mg (see 2) 


Fe = fm(1—mg) (see 1) 
Fs + Fa =} fm (1 — mg) 

Fs + Fa contain each 

; Fs and Fa atoms respectively. 


Q = si— (50 + 2 alk — al) 
This value of Q is obtained as follows: 


Q is the balance of atoms left over after form- 


ing the other bases. 


Total available atoms 
= si+ 100 + al + alk + ti(see Table IX) 








(6) 


(7) 


(8) 


(9) 


(10) 


(a) 
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Total atoms distributed to basis molecules 
3 3 3 3 , 
= 3 alk + 5al+5¢+ 5 alk + 5 fm + ti 


= 3alk+tit 3(alte+ alk + fm) 


Since al + c+ alk + fm = 100 
= 3alk + tit 3 x 100 
= 3 alk + ti+ 150 (b) 
Deducting (b) from (a) 
si —50 — 2 alk + al 
*, Q=si—(50 + 2 alk —al) (11) 
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1. INTRODUCTION 


THE flux and charge spectrum of heavy primary cosmic radiation has been 
studied by several workers by exposing nuclear emulsions at ‘high altitudes, 
and geomagnetic latitudes from 30°N. to 55° N. (Bradt, et al., 1950 a; Dainton, 
et al., 1950; Freier, et al., 1951). To obtain the integral cnergy spectrum 
of the primaries, it is necessary to extend these observations to other lati- 
tudes. The present paper reports experiments made for this purpose near 
the geomagnetic equator. Flights were made from Tambaram, near Madras, 
and from Bangalore, during October and November, 1950. Bcth places 
have geomagnetic latitude 3° N. A summary of the flight data is given in 
Table I. 
TABLE [ 


Summary of Flight Data 





Maximum 
height reached | Plates flown 
in ft. 


Launched | 
trom 


Date of 
launching 


Flight 
number | 


Remarks 





~ 


om O-3192 ore we 


Tambaram 


| Bangalore 
2 


o> 


” 


15-10—1950 | 


20-10—1950 


22-10—1950 | 


29-10—1950 
31-19—1950 
5-11—1950 
8-11—1950 
18-11— 1950 
23-11—1950 


26-11—1950 | 
7-12—1950 | 


8-12—1950 
9-12—1950 


40,000 


76,600 
79,400 
81,300 
50,009 
52,700 
79,600 
81,100 


94,500 
84,700 
91,800 


62,700 





NT4-250u 


G5-500 x» 
NT4-250 uw 
NT4-250 » 

G5-500 uw 
NT4-250 uw 
NT4-250 

G5-500 u 

G5-500 

G5~500 u 


Brass Stack 
Brass Stack 


Brass Stack 


Main line gave way about 50 min. 
aft er launching 


| Plates developed at Bristol 


Plates developed at Bombay 


Plates flown again in’ flight 9 


| Plates flown again in flight 7 
| Plates developed at Rochester 
| Plates lost at sea 


| Highest fight in series, Plates 


lost at sea 


| Plates developed at Rochester 


Balloon break-up. Plates used 
again in flight 13 
Plates not recovered 





A4 
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One object of these experiments was to study the E.-W. asymmetry of 
the heavy primary radiation, which has not yet been measured. For this 
purpose a servomechanism was designed to maintain the plates during flight 
in a fixed orientation. Kodak NT4—-250p and Ilford G5-500y were the 
types used. In a typical flight a stack of 80 plates was used, each 3” by 8”, 
placed vertically with the plane of the plates perpendicular to the magnetic 
meridian. To study the spectrum of the high-energy primaries the “* brass 
stack” was employed. This consists of 20 horizontally stacked plates, 
individually sandwiched between brass plates of 3mm. thickness. The 
plates are positioned with great accuracy. The heavy primaries interact 
principally with the brass, forming narrow showers of break-up products 
which can be followed through successive plates. Either the opening angle 
of such showers or the relative scattering of evaporation alpha particles and 


heavier fragments can be used to determine the energy of the primary (Kaplon, 
et al., 1952). 



























Plates exposed at high altitudes record large 7-meson showers of high 
energy similar to the “‘ R ” shower (Bradt, et al., 19505). The study of the 
development of soft cores in such showers enables one to place an upper 
limit on the lifetime of the neutral 7-meson (Kaplon, et al., 1952). 


The objects of these experiments have not been fully realised, but values 
of the flux have been obtained for the C-N-O group of primaries and for 
the group with Z> 10. The detailed altitude curves have also yielded new 
information on the behaviour of high-altitude balloons. 


2. THE ORIENTATION UNIT 





Fig. 1 shows the orientation unit with the plate bex in position. The 
dome on the right goes over the unit, and is bolted to the circular base. A 
rubber gasket makes the assembly airtight, so that normal pressure is main- 


tained inside during the flight. The dome is of aluminium and is painted 
black. | 








The photograph shows the batteries (A), the electronic circuit (B), the 
plate box (C), the compass box (D), and the monitor (E). Fig. 2 shows 
the circuit diagram. The ends of a compass needle are provided with wires 
dipping into semicircular troughs of water which form part of an impedance 
bridge, to which the output of a 1,000-cycle oscillator is applied. A devia- 
tion of the needle results in a signal whose phase depends on the direction 
of the deviation. The signal actuates a servomotor which rotates the plate 
assembly in such a direction as to annul the deviation. The initial balance 
is brought about by adjusting the two potentiometers. Thereafter the 
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Fig. 2 


orientation of the plate assembly is maintained, provided the rate of rota- 
tion of the gondola is not more than about three revolutions per minute. 
The apparatus showed no tendency to ‘hunt’ or to maintain the gondola 
in oscillation, and with careful preliminary adjustment it behaved well in 
practice. 


The monitor consists of a magnetic needle horizontally pivoted and 
centrally mounted in a flat circular box. The needle carries a polonium 
source at one end, and a circular piece of Kadak NTB 3-35, film is mounted 
inside the lid of the box to record alpha particles from the source. On 
descent the rapid downward acceleration causes the lid to fly open. It is 
then held by a catch, so that no further irradiation of the film occurs to con- 
fuse the record. The distribution of alpha particle tracks on the ,film indi- 
cates the degree to which orientation has been maintained during the flight. 
In flights where the apparatus was recovered intact the monitor record 
showed that orientation had been maintained within about + 15°. A magne- 
tic orientation unit of similar design has recently been described in detail 
(Golian and Reilly, 1951). 
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3. BALLOON TECHNIQUE 


Dewey and Almy balloons were used for these flights, the types being 
K 2,000, J 2,000, J 800 and J 1,400. K balloons are of latex (natural rubber), 
and J balloons of neoprene (plastic). The number indicates the nominal 
weight of the balloon, but there are large individual variations. K_ balloons 
require no pretreatment, but J balloons, in accordance with manufacturer’s 
instructions, are heated for a few minutes in nearly boiling water shortly 
before filling. Many balloons are uneven, and develop bulges or blisters 
on filling, but there is no evidence that these defects appreciably weaken 
the balloons. 


Balloons are filled from hydrogen cylinders through a flexible hose, 
until they just pick up a pre-determined load. Considerable errors are possi- 
ble in judging the precise moment when the load is balanced, particularly 
when there is any wind. Cooling takes place as the hydrogen expands through 
the nozzle, and as the balloon warms up there is a noticeable increase of 
lift. It is therefore difficult to achieve accurately the desired lift, the 
uncertainty being of the order of 10% either way. 


Each flight is planned for a definite rate of rise, given by the empirical 
formula 


Ri see 8 


VE TB 

where R is the rate in 1,000 ft./min., F is the weight the balloon will just 
pick up at the ground, g the external load per balloon during flight, and 
B the balloon weight, all in kg. The formula holds reasonably well for 
balloons attached separately to a long line; but it is not applicable to balloons 
in large clusters, where the rate of rise is noticeably less than given by the 
formula. High rates of rise (800 to 1,200 ft. per min.) appear to be associated 
with a high probability of balloon failure, especially between 35,000 and 
50,000 ft., where turbulent motion is often very evident. 


The necessary hydrogen was produced as required in standard genera- 
tors by means of ferrosilicon and caustic soda. By connecting cylinders 
to the generators it was possible to use double the normal charge. The 
first generator was used to fill one empty cylinder to a pressure of 900 lb./ 
sq.in., and another to 3001b./sq.in. The next charging filled the first 
cylinder to the full pressure of 1,4001b., and the second to 8001b. The 
third charge filled the second cylinder to the full pressure and an empty 
cylinder to 300lb. By repeating this process, five workers could fill 
20 cylinders in about 16 hours, using 15 generators. Previously the generators 
have been used to fill balloons directly, but the present method is more 
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advantageous and involves fewer generators. It does not appear to have 
been used before in India. 


Balloons on 35-foot strings were attached individually or in small clusters 
to a long nylon line at 35-foot intervals. Clusters of 3 to 7 balloons were 
used in different flights. In flight 12, 12 balloons on 80-foot strings and 
11 on 60-foot strings were flown as a single two-layer cluster. Most of these 
survived to a height of 92,000 ft., when atleast a dozen burst simultaneously. 
In this case the central balloons would, when fully expanded, be pressed 
tightly together, and the failure of a single balloon would endanger all the 
others. 


The equipment was attached to the lower end of the line inthe order 
(1) clock release, set to detach the load at a predetermined time, (2) para- 
chute, (3) the gondola carrying the orientation unit, (4) barcgraph, (5) ballast. 
The clock release was used only in the earlier flights. The purpose of the 
ballast is to reduce the initial rate of rise, and to compensate for loss of lif; 
as balloons burst. In the earlier flights 8 kg. of sand, running at the rate 


of 1-2kg. per hour, was used as ballast. In later flights hot water was 
used successfully. 


The launching of a 400-foot line, carrying upto 32 balloons, and with 
equipment weighing upto 40kg., presented some problems. The method 
adopted was to hold both ends of the line, with the top end to windward 
of the load, and allow the central ballcons to lift the line into a vertical arc. 
The top end was then released. When the line was approximately vertical 
the load was released. In this way it was possible to avoid the danger of 
the load being dragged along the ground by an inclined balloon line. 


4. ALTITUDE MEASUREMENTS 


In flight 11 a recording aneroid barometer was used. This gave an 
excellent record, reproducing the details of the trigonometrically deter- 
mined path with great fidelity, but the absolute calibration was uncertain. 
On other flights a specially designed mercury barograph was used, in which 
the movement of the mercury meniscus is photographed on a revolving 
drum. Illumination was provided either by a 2-volt lamp working on vatteries, 
or by daylight. Unfortunately the barograph was tampered with by the 
finders on several occasions, and developed defects on others, with the 
result that only one satisfactory record was obtained, for flight 4, which 
again agrees very well with the trigonometrical trajectory. 


The altitude curves for these flights are therefore based entirely on the 
trigonometrical data obtained from theodolite stations, For the Tambaram 
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flights five permanent stations were used: Tambaram, Meenambakkam 
(6 miles N.E.), Vellore (67 miles W.), Kolar (125 miles W.), and Bangalore 
(177 miles W.). The locations were known with precision, and the three 
distant stations were in continuous touch with Tambaram by two-way short- 
wave radio. Synchronised measurements of altitude and azimuth were 
made at five-minute intervals during the flights. For the Bangalore flights, 
the theodolites at Bangalore and Kolar were used. On the last flight the 
balloons were followed across country by a mobile party taking theodolite 
readings from various locations, making azimuth corrections by reference 
to map positions and sun observations. The resultant altitude curves have 


been computed by an analysis of all the observations available. These curves 
are shown in Fig. 3. 


5. FLOATING LEVELS OF BALLOON SYSTEMS 


The detailed altitude-time curves obtained for these flights have revealed 
an important feature of balloon behaviour which has not been previously 
reported. As the flight rises balloons progressively burst until approximate 
equilibrium is reached. The system then shows a series of stable floating 
levels, which may be maintained for considerable periods. In Fig. 3, flight 
13 shows a single level maintained with very small variations for 9 hours. 


Flights 3 and 4 show cases where the system of balloons falls suddenly from 
one level to another, three such levels being clearly shown in flight 4. When 
ballast is steadily running a ‘level’ becomes a steady increase of height 
with time, but the discrete set of ‘ levels’ is still apparent, as in flights 8 and 
10. The sudden fall from one level to another is almost certainly caused 
by the bursting of a balloon, but owing to the great distances at which the 
balloons were observed this could not always be checked observationally, 


In previous studies it has always been assumed that the pressures inside 
and outside a balloon are equal, and that “the effects of pressure created 
by the elastic contraction of the rubber balloon fabric. ...may be neglected”* 
(Clarke and Korff, 1941). On this assumption the lift due to the hydrogen 
is independent of the pressure. If the total lift is greater than the load the 
balloon rises indefinitely, otherwise it descends to the ground. A group 
would behave in the same way as a single balloon, and could only float if 
the lift happened exactly to equal the load. Even so there could be no float- 
ing level, for there would be equilibrium at all heights. 


The existence of stable floating levels contradicts the above assumption. 

We have to take into account a contribution p, to the internal pressure 
given by 

Py = 2T/r, (1) 
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where T is the tension in the rubber and r the radius of the balloon. At 
the ground p, is of the order of a millibar and is negligible, but it becomes 
increasingly important, both absolutely and relatively, as the balloon expands. 
Rubber balloons blown up on the ground show a minimum value of p, at 
moderate inflation, followed by an increase which is roughly proportional 
to the volume. This suggests that at high altitudes the elastic contribution 
to the internal pressure might reasonably be represented by the equation 

Pr = Py + kip, (2) 


where p is the external pressure and p,, k are constants for a given type of 
balloon. 


We adopt the following notation: 


m = mass of hydrogen in a balloon, 
V=volume of balloon, 
T =the absolute temperature of the surrounding air, 
t = excess internal temperature of the balloon, 
M.y,, M, = molecular weights of hydrogen and air, 
W= total load (= hydrogen + rubber + payload), 
L = W/m = load per gram of hydrogen, 
R = gas constant per mole. 
The volume is given by V = (m/My)R(T + 0/(p+ pi) 
The density of the air is M, p/RT 
The upthrust is thus m(M,/My) (1 + ¢/T)/(1 + pi/p) 
At a floating level this must be equated to W. Hence 
(Ma/My) (1 + t/T)/L = 1 + pi/p 
The value of ¢ is unknown but must be appreciable, since the balloon is 
exposed to the sun’s radiation, and a system of balloons always loses lift 


very quickly at sunset. A fair estimate appears to be ¢/T=0-11. With 
M,/My = 14°4 the equation becomes 


16/L = 1+ p,/p. (3) 
If we assume the validity of equation (2) the equilibrium equation becomes 
16/L = 1 + polp + kip? (4) 


This equation is plotted in Fig. 4 for the values pp=2 mb. and k = 100. 
These values have been chosen tc give a reasonable fit with the observations, 
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Initially a system of balloons has a definite load per gram of hydrogen, 
say L,, which remains constant as long as the balloons remain intact. The 
rise of the balloons is therefore represented in Fig. 4 by the line PQ. AB 
represents the bursting level, at which a balloon reaches its maximum dis- 
tension. At Q one balloon bursts. m grams of hydrogen are lost, and 
W decreases by the weight of the hydrogen and the rubber which falls off. 
L changes discontinuously to L2, and the representative point on the diagram 
is now R. As R is still below the equilibrium level another balloon will 
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burst. The new value of L is Ly, for which the system is above the equi- 
librium level. It therefore falls to S, and if another balloon bursts, to‘T. As 
S and T are below the bursting level the system may persist for long periods 
at these levels. 


The points Q, R, S, T in Fig. 4 have been computed for a flight of 10 
balloons weighing 1-4 kg. each, inflated to a diameter of 1-9 m. or to a lift 
F =2-5kg. For this inflation m=0-3kg. A payload of 16kg., together 
with 10 balloons and the hydrogen, gives a total load of 33 kg. carried by 
3kg. of hydrogen. When a balloon bursts we assume that half the rubber 
(= 0-7kg.) is blown off. The total decrease of W is thus | kg. The burst- 
ing level is assumed to be at 23 mb. With these data we construct Table II. 

TABLE II 








Number of Balloons oat 7 6 
Total load: W a | « 30 29 
Mass of Hydrogen : m xi 0 | 2:1 18 


Load/ kg. of Hydrogen: Lo eel 0 | 14+3 16+1 








Equilibrium pressure in millibars | na a 25-0 | 38-0 | Ground 





The system has therefore two floating levels, and if the bursting level 
were higher there might be 3 or 4. In practice balloons vary in weight, 
some are weaker than others, and it is not possible to fill them all precisely 
to the same lift. The amount of rubber lost on bursting is also very vari- 
able, and there are cases where balloons develop pinholes and shrink slowly, 
so that no rubber is lost at all. For these reasons floating levels cannot be 
predicted accurately in advance, but the present results show that the flights 
behave substantially in the manner to be expected from this theory. 


The rise in p, for high inflations does not occur with neoprene balloons 
when blown up to bursting point at the ground though the rubber balloons 
show it clearly. The difference is almost certainly a temperature effect. 
The manufacturers, in a private communication, have stated that neoprene 
‘freezes’ or becomes less estensible at stratosphere temperatures. This 
would correspond to a large value of k in the above theory. 


6. DETERMINATION GF PRIMARY FLUX 


The experiments provided 4 sets of plates for analysis, obtained in 
flights 3,4. 7 and 12. The plates of the brass stack in flight 12 showed strip- 
ping of the emulsion on development, and have not yielded satisfactory 





High Altitude Balloon Experiments 51 


measurements. The mean altitude in flight 3 was rather low, and though 
the plates are good, the primary tracks are too few to give a reliable deter- 
mination. The results reported are therefore based on flight 4, and flight 7. 
Flight 4 was analysed at Bombay and flight 7 at Rochester. 


Of the 80 plates from flight 4, the two end plates and two from the middle 
were scanned systematically for heavy tracks. When found, these were 
followed from plate to plate through the stack to distinguish primaries from 
slow secondaries. Primaries have an energy above the cut-off, which is 
approximately 6-5 Bev/nucleon at the geomagnetic equator, and at these 
energies the primaries pass completely through the stack unless they suffer 
nuclear collisions. The mean free path in the material of the stack is known 
to be 34-5 g./cem.* for incident particles with 6< Z< 10, and 25-5 g./cm.* 
forZ> 10. A break-up occasionally occurs in the emulsion, and the break- 
up products can often be recognised as alpha particles (4 times minimum 
ionisation) and singly charged particles (minimum ionisation). In such 
cases, a lower limit to the charge of the primary may be assigned from the 
break-up products. 


5-Ray measurements were made on all tracks for which Z could be 
assigned in this way with some assurance, and also on alpha particle tracks. 
The number of 5-rays in a standard length of track (100) is given by the 


formula (Bradt and Peters, 1950 c). 
Ns = AZ? + B. 


From the results obtained from the known tracks, the constants A and 
B were evaluated, and the relation was then used to determine the charge 
on all the primary tracks found. As all the primaries have relativistic speeds 
the 5-ray count over a reasonable length of track, usually 1,000 or more, 
gives the value of Z. The uncertainty is cf the order of one unit in Z. For 
analysis the tracks are divided intc two groups (1) those for which 6< Z< 10 
(generally called the C-N-O group), (2) those for which Z> 10. For each 
track the orientation and dip (with due correction for emulsion shrinkage) 
were recorded, these angles being then converted to true zenith angle and 
azimuth. The direction of travel is determined from the fact that a primary 
particle must be moving downwards. 


Let z be the vertical axis normal tc the plate. The zenith angle of the 
track is 6, and its azimuth, measured from the x,z plane, is 6. The dip 8 
is the complement of the angle between the track direction and x. These 
angles are related by the equation 


sin B = sin 4 cos ¢. 
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Owing to difficulties of recognition, only tracks of a certain minimum 
projected length are included in the analysis. This is equivalent to the 
condition B < By. This implies that 4 > 4,4, 


where 
cos $9 = sin B,/sin 0. 


When @ is less than By, all values of ¢ between 0 and z/2 are included, but 
when @ exceeds f, the permissible azimuth ¢ has the lower limit 5. 


Consider an isotropic incident flux of N, particles per sq. cm. per sec. 
per sterad. The plate presents a projected area (A sin f), where A is the 
area of the plate, and the element of solid angle is sin 6d@d¢. The rate 
of incidence on the plate at the angles 6, ¢ is therefore 


n(6, 6) = N, Asin? 6 cos ¢ dé dd 


Integrating over ¢ with the given limits, the rate of incidence for a given 
zenith angle becomes 

n(6)=2N,A/S (8) dd, 
where 

f (9) = sin?0 for 6< By 


f() = sin? 0(1—,/1 — S28 = nats) for 8 > By. 


To allow for the absorption in the atmosphere, n (6) must be multiplied 


by a factor e~"/*°°S® where A is the mean free path in air, and h the vertical 


height of the atmosphere above the point of observation. h is known as 
a function of the time from the altitude curve, and the total number of tracks 
recorded at the zenith angle @ is thus given by 


T —h/X cos 0 
No = 2N,Af(6) d0fe dt. 
ry 


Denoting the integral by ¢ (6), the total number of tracks recorded between 
any given limits of zenith angle is 


6 
N (64:5 65) = 2NoA ff (0) $ (6) a8. 


Some of the observed tracks represent particles which have had an 
appreciable path in the stack before reaching the selected plate. If the path 
in the stack is L, and the mean free path A’, cach track must be assigned a 
statistical weight el’. The appropriate value of N (4, 65) to employ is 
the sum of the statistical weights of all the observed tracks. 


The analysis has been made for several plates of the stack, and for 
different values of By. High values of By give smaller flux values, showing 
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that tracks (especially of the C-N-O group) are liable to be missed or excluded 
by the observer if the minimum projected length is too small. The value 
of By finally adopted is 34°, corresponding to a minimum projected length 
of 300u in a 200p emulsion. With this criterion a single end plate has 
yielded 24 tracks, with a total statistical weight of 33. 


Taking the mean free path in air of the C-N-O group as 27 g./cm.? (Bradt 
and Peters, 1950c), ¢(9) is determined by treating the altitude curve as a 
series of sections of constant height. The product f(6) 4 (@) is then tabulated 
at 10° intervals, and the integral is evaluated by quadrature. The procedure 
for the group Z > 10 is similar, except that the value 21 g./cm.? is used for 


the mean free path in air, and the value 25-5 g./cm.? for the mean free path 
in the material of the stack. 


The tracks recorded in flight 4 show an apparent excess of particles 
coming from the south. This is almost certainly due to an unforseen tilt 
of the stack, which would cause an excess of particles to be recorded on the 
face tilted upwards. A tilt of some 10° would account for the observed 
asymmetry. Calculation shows that the total number of particles recorded 
in the tilted plate will not differ, within the limits of accuracy of this experi- 
ment, from the number recorded in a vertical plate. 


The final values adopted, after considering these results and compar- 


able results of flight 7, are as follows: 
For 6<Z< 10, Ny = 1-45 + 0-30 particles/m.? sec. sterad 
For Z> 10, N, = 0:33 + 0-08 particles/m.? sec. sterad. 


7. CHARGE AND ANGULAR DISTRIBUTION 


The charge specturm of the primaries has been determined by others 
at higher latitudes (Bradt and Peters, 1950a; 1950c). The present results 
agree with this in general, but some uncertainty in the 5-ray calibration and 
the relatively small number of particles observed make the specturm un- 
certain in detail. The relative intensity of the two groups of particles agrees 
well with previous observations. 


For a major part of the duration of the flight, the plates remained under 
about 55 g./cm.? of matter. The primaries have large collision cross-sections, 
and hence those at large zenith angles were completely cut off. This resulted 
in a relatively small observable E.W. asymmetry of the incoming heavy pri- 
mary radiation. Thus, this asymmetry, in the present experiment, is within 
the statistical uncertainty of the results. 
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8. SUMMARY 


The paper describes high-altitude balloon experiments with nuclear 
plates, made near the geomagnetic equator. A magnetic orientating unit 
was used. The behaviour of balloons is discussed, and it is shown that 
systems of balloons may have a series of stable floating levels at determinate 
altitudes. Analysis of the tracks observed in the plates has given the follow- 
ing values for the primary cosmic ray flux at the geomagnetic equator: 


For 6< Z<10, N, = 1-45 + 0-30 particles/m.? sec. sterad. 
For Z> 10, N, = 0:33 + 0-08 particles/m.? sec. sterad. 


Any observable E.W. asymmetry at the heights reached by the plates is 
within the statistical uncertainty of the results. The charge spectrum is 
in general similar to what has been previously observed. 
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1. INTRODUCTION 


In Part I, the author published the results of a study on the infra-red absorp- 
tion by a single crystal of potassium chlorate in the range 1 to 15y, by 
means of a Beckmann infra-red spectrophotometer (model I R2). In the 
present investigation the work has been extended upto 22,» using thinner 
crystals than formerly. The instrument itself had meanwhile been over- 
hauled and its sensitiveness improved. The results now reported are there- 
fore superior to those presented in Part I. 


The crystals were studied, as before, with the flat face (001) perpendi- 
cular to the infra-red radiation, percentage transmission being read off 
directly with the instrument. Readings were taken at 0-1, intervals and 
at still closer intervals near absorption maxima. The percentage cut-off 
curves of two of them which show the absorption maxima most clearly 
have been reproduced in Fig. 1. Curve III represents the % cut-off of the 
specimen of thickness 0:15mm. and Curve I of specimen of thickness 
0-06 mm. 

The following features are shown by the curves: 


I. A sharply defined, strong absorption band with its maximum at 
5:32 (1887 cm.—). 

Il. Two well-defined absorption maxima at 6-3 (1587 cm.) and 
6-9 (1449 cm.) which are weak in comparison with the 5-3 band. In 
Part I of the paper, this doublet was not resolved and a rather broad band 
with its maximum at 6-7 was recorded. Coblentz (1906) to whom no 


teference was made in the previous paper through oversight had noticed 
these bands as also the one at 5-3. 


Ill. A very strong absorption band with a flat top, extending from 
10n to 10-8 (1000cm.-' to 926 cm.-?). 


IV. A medium strong absorption band at 13-5, (741 cm). 
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V. A strong and sharply defined absorption maximum at 16-2, 
(617 cm.~). 


VI. In the slope of the absorption maximum at 16-2, we notice 
inflexions at 17-25 (580 cm.—), 18 (556cm.~*), and 18-75 » (533 cm), 


Vil. A strong and sharply defined absorption maximum at 20-3, 
(493 cm.-'). This band shows indications of structure. 


Schefer and Schubert (1921) in their reststrahlen studies with pressed 
powder of KCIO, found reflexion maxima at 10-12 and 16-12, which 
evidently are related to the absorptions referred to in III and V above. 
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Fic. 1. Infra-red absorption spectra of single crystals of potassium chlorate. 


2. INFRA-RED ACTIVE FUNDAMENTALS 


In Part I, the absorption maxima in the infra-red were explained in 
a general way in terms of the characteristic frequencies of the chlorate ion 
the free state, viz., v,(930cm.-*) and », (620 cm.-) belonging to the totally 
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symmetric class A,; vg (975cm.') and v, (486cm.-') belonging to the 
doubly degenerate class E. But in the transformation from a free ion to a 
.crystal, the number of internal vibrations increases because there are two 
molecules in each unit cell and the symmetry of the crystal structure 
(C,,”) removes the degeneracy of modes of the E class. The group analysis 
gives the number of internal vibrations and the selection rules as shown, 
in Table I. 
TABLE | 


Number Raman effect | Infra-red 


Active Inactive 
Active Inactive 
Inactive Active 


Inactive Active 


ie., there are six Raman active internal vibrations and six infra-red active, 
the two sets being mutually exclusive because of the presence of a centre 
of symmetry. We can picture the change in internal frequencies from 
a free ion to crystal as follows: 


TABLE II 








Free ion Crystal 
Raman active 
and Infra-red Raman active 
active 





v, (A,) v, (A;) | w, (Bs) | 
vy (Ay) v3 (A;) | ws (Ba) 

vy (EB) vy (A,), v2’ (As) | w, (Bs), w,’ (B,) 
v,(E) vy (Ay), v4 (Az) | w,(B,), 4’ (B,) 


| 
| 
} 
| 





Thus the infra-red and Raman frequencies would, strictly spreaking, be 
different. 


AS 
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The Raman effect in the powdered crystals has been studied by Krishna- 
murthi (1930) and Venkateswaran (1938), and by Shantakumari (1950) in 
a single crystal. Their results have been given in Table III. 


TABLE Ill 








Krishnamurthi | Venkateswaran Shantakumari C10.- ion 
(Powder) (Powder) (Single crystal) | . 





915 (1) | 920 (1) 
929 (10) 930 (10) 
977 (4) | 975 (7) 
615 (1) | 620 (1) 





478 (3) | 486 (6) 
493 (0) 


* With two faint companions. 


We notice that the Raman frequencies which have been reported fall 
into four main groups which approximate respectively to the four funda- 
mental frequencies of the ClO, ion. The results of the different authors 
are somewhat divergent, which may in part be ascribable to their having 
recorded the spectra in different circumstances. For our present purpose, 
we may ignore the splitting of the fundamental frequencies in the crystal 
shown in Table II and base our interpretation of the infra-red absorption 
data on four fundamental frequencies denoted as w,, wa, w3, wa respec- 
tively, for which we assume the values 930, 975, 620 and 486 cm.~!, these 
being the same as those observed in the Raman effect by Shantakumari. 


3. INTERPRETATION OF RESULTS 


On the assumptions made above, the following interpretations of the 
observed bands are put forward as the most probable (Table IV). 


The bands IV and VI have been interpreted as combinations of the 
internal and lattice frequencies of the crystal. There are actually nine 
Raman active and six infra-red active lattice frequencies, the two sets being 
mutually exclusive. The selection rules permit combinations between infra- 
red active internal frequencies and Raman active lattice frequencies and 
vice versa. The infra-red active lattice oscillations are unknown. Shanta- 
kumari has recorded five lattice lines in Raman effect which we may denote 





Infra-Red Absorption Spectrum of Potassium Chlorate Crystals—II 59 : 
TABLE IV 





Observed Absorption Maxima Calculated 
Assignment 


Au cn:* Au cm.~? 














I 5°3 1887 5-25 1905 


II 6-3 1587 6°27 1595 
6-9 1449 6°84 | 1461 


Il 10-10-8 1000-926 10-26 975 
10-75 | 930 


IV 13-5 741 13-4 747 
Vv 16-2 617 16-13 | 620 


VI 18-0 556 17-7 566 
18-75 533 18-5 540 


Vil 20-3 493 20-58 486 


























as L, (54 cm.-*), L, (82 cm.-), L, (98 cm.-'), L, (127 cm.) and L; (145 cm). 
Since the data regarding the lattice oscillations are incomplete, the assign- 
ments shown against IV and VI are necessarily somewhat conjectural. 


Some interesting features are noticed as regards the relative intensities 
of the internal vibrations in Raman effect and infra-red. In Raman effect, 
all the workers have got 620 cm.~! line.very weak. In fact, Krishnamurthi 
in his study of Raman effect in KCIO, crystal powder and Venkateswaran 
in Raman effect in solution of KCIO, failed to observe this line. On the 
other hand, corresponding to that line we notice a strong and sharp absorp- 
tion in infra-red. The 486 cm.— line which is strong in Raman effect is also 
strong in infra-red absorption. As regards 930cm.-! and 975cm.-', the 
indications are that the latter is more strongly active than the former in 
absorption, while in Raman effect, 930cm.—' is stronger than 975 cm.-}, 
though it has a comparable intensity. 


Iam very grateful to Prof. Sir C. V. Raman, F.R.S., N.L., for the encourage- 
ment and guidance he gave me during this investigation. My thanks are 
also due to the Investigator-in-charge C.S.I.R. Radiation Scheme, Pgona, 
for his kind permission to use the infra-red spectrophotometer. 
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4. SUMMARY 


Crystals of potassium chlorate exhibit strong absorption bands in the 
infra-red at 10y-10-8y, and at 16:2 and 20-3. These are identified 
with the active fundamentals at w, (930 cm.-!, 10-75 pn), w, (975 cm, 10-26 p), 
ws (620cm-1, 16-134) and w, (486 cm.-!, 20-58). It is noticed that w, 
is quite strong, while the corresponding Raman frequency comes out very 
weakly. The vibration w, is more strongly Raman-active than w., while 
the reverse appears to be the case in the infra-red. Absorptions are also 
observed at 5-3, 6-3 and 6-94 which are interpreted as first-order 
combinations among these four fundamental frequencies, while other bands 


at 13-5, 18 and 18-75 appear to be combinations of internal frequen- 
cies with low frequency lattice vibrations. 
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INTRODUCTION 


It is a common experience in the sea in low latitudes that insolation warms 
the upper layers, and that below a depth of about 100 meters there is usually 
a rapid fallin temperature. This lower layer of water with rapidly increasing 
density is known as the thermocline or discontinuity layer. The presence 
or otherwise of a thermocline has far reaching effects on the organic produc- 
tion and the distribution of planktonic organisms. 


The earliest account of the temperatures at different levels of the Indian 


Ocean was given by Carpenter (1887) in which he gives the results of his 
observations in the Bay of Bengal. Schott (1902), Kriimmel (1911), Matthews 
(1926) and Sewell (1932) have added to our understanding of the temperature 
gradients in the waters surrounding India. Still our present day knowledge 
on this subject is meagre. 


Many devices have been invented for obtaining the temperature as a 
continuous function of depth but the bathythermograph designed by Spilhaus 
(1938) and later considerably modified by him, has overcome many of the 
difficulties that rendered previous designs ineffective. With the bathy- 
thermograph (Thermarine Recorder Model FA-190011) available at the 
Central Marine Fisheries Research Station it was decided to make observa- 
tions on the temperature gradients whenever possible. Although the aim 
is to obtain a complete picture of the kaleidoscopic changes in the hydro- 
logical conditions throughout the year, owing to many difficulties it wil] 
only be possible at this stage to make intermittent observations and try to 
complete the picture in course of time. 


* Published with the permission of the Chief Research Officer, Central Marine Fisheries 
Research Station, Mandapam Camp, South India. 
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The following observations were made during a cruise from Trincomalj 
to Madras, Calcutta, back to Madras and then to Mandapam. During the 
cruise observations on temperature gradients were made at 22 stations and 
water samples were collected at these stations as well as from another five 
stations (Fig. 1). The pH of these samples was determined on board the 
ship using Bromthymol blue as the indicator but the determination of salinity 
had to wait until the samples were brought back to the laboratory. The 
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Fic. 1. Shows the stations at which observations were made. Numbers underlined 
indicate stations where bathythermograph readings were not taken. 
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salinity was estimated by the standard chlorinity titration method. Informa- 
tion such as the exact position, atmospheric temperature, barometric pressure, 
golour of water, nature of sea, depth and nature of bottom, direction and 
force of wind at the time of observation was also recorded in log sheets. At 


stations covered during day time, the penetration of sunlight was measured 
with the aid of a Secchi disc. 


It may be mentioned here that the observations presented below are 
not to be regarded in any sense as a complete survey, but represent the first 
step in collecting adequate data for an understanding of the hydrological 
conditions of the water masses surrounding our peninsula. 


In embarking on investigations of this kind it is necessary to call upon 
many persons for assistance and it gives me great pleasure to acknowledge 
the help received. This investigation, carried out on board I.N.S. Delhi, 
was made possible by the willing co-cperation of the Navy and the Central 
Marine Fisheries Research Station. First and foremost amongst these, 
I wish to mention Rear-Admiral G. Barnard, Rear-Admiral Commanding 
the Indian Naval Squadron and Captain Chatterji to whom I tender my 
gratitude and appreciation. To the Naval Personnel who helped me in 
managing the bathythermograph and the Secchi disc I offer my sincere thanks. 


DISTRIBUTION OF TEMPERATURE IN THE UPPER 200 FEET OF WATER 


Carpenter (1887) showed that from the surface to a depth of nearly 
20 fathoms the temperature shows only little change but below this lies a 
zone extending downwards for more than the next hundred fathoms in which 
the temperature falls rapidly. Schott (1902) concluded from the data he 
had at his disposal that in the Indian Ocean the discontinuity layer lies at 
a depth below the surface which is intermediate between that found in the 
Atlantic and Pacific Oceans, i.e., between 90 and 140 meters. He attributes 
the presence of this zone partly to the alternating effects of the seasons and 
partly also to the effect of lateral movements of the upper stratum which 
cause a thinning of the surface currents resulting in the deep lying colder and 
more saline water approaching the surface. Thus during hot weather the 
surface water becomes condensed and has a raised temperature but during 
the cold season this water cools and sinks down to a depth of 91 to 137 
meters. Kriimmel (1911), however, points out that true discontinuity layers 
in the other seas occur at a level that is only a few meters, 15 to 20 meters, 
below the surface which is caused by the heating up of the uppermost levels 
in summer and refuses to believe that this sudden change in temperature 
at depths as great as 100 meters can in any way be due to the heating up and 
¢vaporation of the surface water.’ He is of opinion that a sharply defjned 
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thermocline is caused either due to the ascending movements of deep strata 
that carry cold water towards the upper levels, which have been warmed by 
the sun, or to horizontal movement of large masses of water due to the surface 
currents, the warm upper layer thus becoming sharply separated’ from the 
deeper colder water. Sewell (1932) suggests that in the Bay of Bengal the 
discontinuity layer may be produced in yet another manner, at least at certain 
seasons of the year, by the dilution of the surface water by rain or river water 
and thus by virtue of its reduced salinity and specific gravity it floats on the 
top of more saline water preventing usual convection currents setting up and 
restricting warming up by insolation tc the upper stratum resulting in the 
formation of a discontinuity zone where two layers of unequal salinity meet. 
According to him the average temperature of the surface water, as calcu- 
lated from all observations is 28-31°C. and, as a rule, there is only very 


little change in the upper 25 fathoms at which depth the average temperature 
is 26-75°C. 


Sewell, op. cit., has constructed a graph showing the upper and lower 
limits of the discontinuity zone throughout the year fer the Bay of Bengal 
and the Andaman Sea. He also points out, however, that no data were 
-available for the months of either August or September. The upper level, 
according to Sewell, for August is 40 to 50 fathoms and the lower limits for 


the same month being from 78 to 90 fathoms. But the observations recorded 
here for the month of August (August 9 to 29) show the existence of definite 
thermoclines at several places in the Bay of Bengal at comparatively low levels. 


Schott (1902) has suggested that the thermocline should be defined as 
one in which the temperature changes by not less than 2°C. (3-6° F.) in 25 
meters vertically. According to this definttion it will be noticed that there 
are well marked thermociines developed within the upper 200 feet at several 
places during August (stations 8 to 11, 17 to 20, 25 and 26). The depths at 
which thermoclines develop vary from place to place and also from time to 
time. It is seen from Figs. 3 and 4 that at the stations 8 to 11 thermoclines 
developed approximately at 65, 65, 70 and 75 feet respectively, whereas at 
stations 17 to 20, 25 and 26 the thermocline was at 130, 90, 35, 60, 120 and 
110 feet respectively (Figs. 8, 9, 11 and 12). It may be added here that at 
stations 12 to 16 (Fig. 1) one would not normally expect definite thermoclines 
in the upper 200 feet of water as there is likely to be considerable mixing up 
of water in the upper layers owing to the influx of large quantities of water 
brought down by the Ganges and Mahanadi (Figs. 5, 6, 7 and 8). In this 
connection it may be pointed out that the depths at which thermoclines deve 
loped at stations 8 to 11 gradually increase from 65 to 75 feet which may 
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Fics. 2-7. Show the distribution of temperature in different depths at ten stations. 
Number against each line represents station number. 


indicate the spreading of a tongue of less saline and lighter water which, as 
is to be expected, thins out gradually. It is likely that thermoclines are deve- 
loped at the stations mentioned above at the places where the two strata of 
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Fig. 10 Fig.l Fig.12 


Fics. 8-12. Distribution of temperature at different depths at twelve stations in the 
Bay of Bengal. Number against each line denotes station number. 


water of different salinity meet, as suggested by Sewell, and it is probable 
that thermoclines might have existed at deeper layers as we approach the 
mouths of the rivers, Again the thermocline at station 17 was at 130 feet, 
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the location of the station being near the mouth of Mahanadi, whereas as 
we move away from the river mouth it will be observed that the thermocline 
develops closer to the surface—90 and 35 feet at stations 18 and 19 res- 
pectively. Station 20 shows a thermocline at 60 feet while at the next four 
stations there are no thermoclines developed within the upper 200 feet 
(Figs. 9, 10 and 11). Thermoclines are developed at 120 and 110 feet 
at stations 25 and 26 respectively. The same explanation holds good for 
all these stations, i.e., those located just opposite river mouths (21 to 24) 
do not show thermoclines within the 200 feet and at those places where the 
influence of river water is only slightly felt the thermoclines develop compara- 
tively near the surface. It is of interest to note here that at stations 8, 10, 
18 and 19 the difference in the temperatures between the surface and at 200 
feet is about the same as the maximum difference observed in the fluctuation 
of the surface temperature in any month. According to Sewell (1929) the 
extreme range of surface temperature, i.e., between the highest and the 
lowest recorded temperatures, in each month rises from 3-2°C. in April 
to 6:-45°C. in December. The range in surface temperature recorded during 
this.cruise is 4-02°C, (ca., 7-23°F.). Stations 8 to 11 and 18 to 20 show 
this or greater difference in temperature between the surface and at 200 feet 
(Figs. 3, 4, 8 and 9). This may mean that the changes in the surface tempe- 
rature may be brought about by the mere upwelling of the subsurface water 
but whether this phenomenon actually happens and also whether this is one 
of the contributing factors it is difficult to say at present. 


LIGHT PENETRATION 


Another factor which is most important for organic production is radia- 
tion. It not only directly controls the phytoplankton production but is 
important indirectly through its effect on temperature. The light that is 
effective for phytoplankton growth depends, of course, on the transparency 
of the water. The light falling on the sea is continuously changing with time 
and thus the depth of photosynthetic zone varies throughout the day with 
the intensity of illumination. The transparency of the water plays a leading 
part in controlling the depth of the photosynthetic zone and similarly the 
relative position of the thermocline and the euphotic zone is also of consi- 
derable importance in organic production. 


Factors such as the presence of debris and plant cells in large quantities 
increase the turbidity. In certain cases plant cells themselves may become 
abundant to cause a significant decrease in transparency. The extent to 
which plant cells may cut off light is shown by Clarke (1946) in which he 
has given the relation between the Secchi disc measurements and the abund- 
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ance of phytoplankton. In all cases in which plant pigments exceeded 
40,000 Harvey units per m* the Secchi disc depth was found to be reduced 
to 8 meters or less. 


In order to get an idea of the penetration of light in the Bay of Bengal 
limits of visibility tests were conducted with the aid of the Secchi disc. It 
will be noticed from Table I that the depth varies from 33 to 78 feet except 
at station 15. This station, as will be seen from Fig. 1 is located very close 
to the mouths of the Ganges and due to the influx of large quantities cf mud- 
laiden water from the river, the transparency of water was reduced. Conse- 
quently the Secchi disc reading was only 14 feet. Owing to lack of facilities 
it was not possible to make plankton collections at any of these stations and 
hence the organic production in the euphotic zone could not be estimated. 
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r has already been reported that the root powder (moisture content: 13-42 

per cent.) of Asparagus racemosus Willd. contains 37-44 per cent. of free 
sugars and 2-10 per cent. of mucilages, while the total non-fibre carbo- 
hydrates amount to 52-89 per cent.1. The other constituents of the non- 
fibre carbohydrates have now been examined and are found to be hemi- 
celluloses (5:40 per cent.) and some insoluble polysaccharides (7-95 per 
cent.). The isolation and examination of the hemicelluloses are reported 
in this paper. 


After the removal of saponins, free sugars, mucilages and pectins by 
successive treatments with ethyl acetate, 70 per cent. alcohol, hot water and 
ammonium oxalate solution, the root powder has been extracted for the 
isolation of the hemicelluloses with 4 per cent. sodium hydroxide according 
to the method of Norman.? On further fractionation as per the procedure 
of Norris and Preece,* the hemicelluloses have given rise to two fractions, 
namely, A, and B,. 


Both A, and B, are amorphous substances, which are almost white in 
colour with a greyish tinge. They are soluble in hot water, the solutions 
being slightly viscous but not mucilaginous. Further, they are readily hydro- 
lysed by boiling dilute mineral acids, indicating their hemicellulosic nature.‘ 
They are composed of the same sugar and uronic acid, namely, glucose and 
galacturonic acid, but in different molecular proportions. In A, the mole- 
cular ratio of glucose to galacturonic acid is 10:1, while in B, it is 5:2. 


The two hemicelluloses have been subjected to partial hydrolysis also. 
This degradation is best brought about by heating the material with 1-5 
per cent. sulphuric acid for 3 hours on a beiling water-bath. It may be 
noted that direct boiling of the mixture under reflux even for a short time 
usually results in almost complete hydrolysis. Under the controlled condi- 
tions mentioned both A, and B, are hydrolysed to glucose and an aldobionic 
79 
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acid. The latter has been isolated as its barium salts On further hydro- 


lysis it gives rise to d-glucose and d-galacturonic acid in equimolecular 
proportions. 


It may be noted that the aldobionic acid present in the mucilage portion 
of these tubers is also a gluco-galacturonic acid.® It is interesting to note 
that in both the mucilage and the hemicellulose portions the aldobionic 
acids are composed of the same sugar and uronic acid; however, they differ 
in their resistance to hydrolysis by acids. The constitutional factors res- 
ponsible for this difference are yet to be studied. 


It is also interesting to note that amongst the three species of Asparagus, 
the tubers of which have been examined in some detail, A. filicinus does nct 
seem to contain any hemicelluloses,® while the other two contain two frac- 
tions each.’ Further, A. adscendens and A. racemosus differ as regards the 
nature of the sugars and also the uronic acids constituting the hemicelluloses. 
The hemicelluloses from A. racemosus are composed of glucose and gala- 
cturonic acid, while those of A. adscendens are constituted from xylose, 
glucose and glucuronic acid. 


EXPERIMENTAL 


Isolation and fractionation of the hemicelluloses—The debarked root 
powder of Asparagus racemosus was extracted successively with ethyl acetate, 
70 per cent. alcohol, warm water and 0-5 per cent. ammonium oxalate solu- 
tion in order to remove saponins, free sugars, mucilages and pectins. The 
residual powder was worked up for the isolation of the hemicelluloses. It 
(100 g.) was heated with 4 per cent. sodium hydroxide (500 c.c.) at 45-50° C. 
for two hours with frequent shaking. The mixture was filtered through 
a fine muslin and the residue was twice again extracted with the alkali, taking 
the same volume each time. After clarification by repeated filtration through 
glass wool, the extract was treated with 10 g. of sodium hypochlorite and 
just acidified with hydrochloric acid. This treatment was to destroy the 
lignins, if present to any extent (method of Norman’). After 10 minutes 
when no more chlorine was evolved and when the hemicelluloses settled 
down, the supernatant liquid was decanted off and the rest centrifuged. 
The separated solid was washed first with small quantities of water, then 
with 60 per cent. alcohol and finally with hot absolute alcohol. It was then 
dried first in air and then in a desiccator. The yield was 5-40 per cent. on 
the weight of the debarked root powder taken. 


The product was fractionated according to the method of Norris and 
Preece.? A 1 per cent. solution of the substance in 4 per cent. sodium 
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hydroxide (500 c.c. in volume) was treated with excess of glacial acetic acid 
and allowed to stand for 6 hours. The solid that settled down (Fraction A) 
was separated in a centrifuge, and washed successively with water and hot 
alcohol. Its yield on the basis of the root powder was 2-00 per cent. The 
mother-liquor (600 c.c.) was next treated with half its volume of ethyl alcohol, 
when Fraction B separated out as a greyish white substance. It was also 
isolated, washed and dried as Fraction A, and was obtained in 3-40 per cent. 
yield. On further treatment with excess of alcohol, the mother-liquor, 
left after the separation of B, did not give any more solid (Fraction C 
absent). 


The above fractions were subjected to further fractionation. They 
were separately dissolved in hot 4 per cent. sodium hydroxide so as to form 
1 per cent. solution and filtered through glass wool to remove the slight turbi- 
dity present. On treatment with Fehling’s solution at the rate of 30c.c. 
for every 100c.c. of the alkaline solution, no precipitate separated out in 
either case, indicating the absence of Fractions A, and B,. Hence the two 
hemicelluloses corresponded to Fractions A, and B, of Norris and Preece. 
On acidification, followed by the addition of excess of alcohol, the alkaline 
solutions precipitated the hemicellulose. Both A, and B, were purified by 
repeated dissolution in dilute alkali and reprecipitation by means of acid 
and alcohol. Three such treatments were required to get ashless products. 
The final yields of purified A, and B, were respectively 1-65 per cent. and 
2:90 per cent. on the weight of the debarked root powder. 


Hemicellulose A,.—In the purified condition A, was greyish white in 
colour, and contained 9-26 per cent. of moisture. It was insoluble in alcohol 
acetone, ether, etc., but was soluble in hot water yielding viscous but non- 
mucilaginous solutions. In 0-5 per cent. aqueous solution its specific rota- 
tion at 18°C. was + 19-7°. 


The hemicellulose contained uronic acid, as it responded to the naphtho- 
resorcin test.* The uronic acid was quantitatively estimated according to 
the method of Dickson, Otterson and Link® and was found to be 8-8 pef 
cent. on Zero-moisture basis of the hemicellulose. It did not contain any 
methoxyl group or pentose. 


The purified material (2 g.) was hydrolysed by boiling under reflux with 
2 per cent. sulphuric acid (150 c.c.) for 2 hours. The hydrolysate was neutra- 
lized with barium carbonate, filtered and the filtrate examined for the sugars 
and the uronic acids adopting filter-paper chromatography (horizontal migra- 
tion method of Rao and Beri.!°) Only glucose and galacturonic acid could 
be detected. The identity of glucose was confirmed by the preparation of 
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its osazone with its characteristic crystal structure and melting point (204- 
206°C.). The identity of the uronic acid was also confirmed by the deter- 
mination of the specific rotation of its barium salt {[a]%°° = + 24-8°} 
and its oxidation with nitric acid of sp. gr. 1-15 to produce mucic acid. 


Composition of Hemicellulose A,.—Since the hemicellulose contained 
only glucose and galacturonic acid, and the latter was found to be 8-8 per 
cent., the composition of the hemicellulose might be taken to be glucose 
and galacturonic acid present in the molecular ratio of 10:1. 


Hemicellulose B, and Its Composition.—This fraction too was similar 
to Fraction A, in its colour, amorphous structure and solubility. It con- 
tained 8-2 per cent. of moisture. Its specific rotation in 0-3 per cent aqueous 
solution was + 27-1° at 18°C. It also underwent easy hydrolysis with 
dilute boiling sulphuric acid. The hydrolysate was examined as in the case 
of hemicellulose A,, and was found to contain the same sugar and uronic 
acid, namely, glucose and galacturonic acid. In this case, the sugar and 
the uronic acid were present in the ratio of 5:2. 


Partial Hydrolysis of A, and B,.—The purified hemicellulose A, (5 g.) 
was heated in a round-bottom flask with 1-5 per cent. sulphuric acid (150 c.c.) 
on a boiling water-bath (98° C. at Dehra Dun altitude), when it rapidly went 
into solution. The course cf the hydrolysis was followed by noting the 
optical rotation at frequent intervals. After 3 hours, when there was no 
more change in the rotation, the solution was neutralised with barium carbo- 
nate in the hot, concentrated to a small volume (40c.c.), filtered and then 


treated with 95 per cent. alcohol (100c.c.). The precipitated solid and the 
filtrate were examined separately. 


From the filtrate alcohol was distilled off under reduced pressure. The 
resulting aqueous solution did not show the presence of any oilgosaccharides, 
since on further hydrolysis it did not undergo any change in reducing power 
or rotation. It contained only glucose. 


The precipitate (barium salt) was purified by repeated dissolution in 
water and precipitation by means. of alcohol. After the third precipitation 
it was filtered, dehydrated using hot absolute alcohol and dried in a vacuum 
desiccator. It contained 16-40 per cent. of barium (estimated as barium 
sulphate) and liberated 9-89 per cent. of carbon dioxide on heating with 12 
per cent. hydrochloric acid. These data indicated that the substance 
was the barium salt of an aldobionic acid. [Barium aldobionate, 
(C,,H,,0,,COO), Ba, contains 16-21 per cent. of barium and liberates 10-38 
per cent. of carbon dioxide]. The salt (2g.) was further hydrolysed by 
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boiling with 3 per cent. sulphuric acid (200c.c.) under reflux for 2 hours) 
As the separated barium sulphate did not interfere with the course of boi: 
ing, no attempt was made to remove it. On analysis the hydrolysate wag) 
found to contain only glucose and galacturonic acid in equimolecular Pros 
portion. s 


B, was also hydrolysed in the same way, and the same products were 


obtained. * 
SUMMARY Be 


The hemicelluloses of the tubers of Asparagus racemosus Willd. have 
been isclated and chemically examined. Two fractions, viz., A, and B, ¢ 
Norris and Preece designation, have been obtained. 


Both the fractions are constituted from the same sugar and uronic acid, 
namely, glucose and galacturonic acid but in different proportions. 
hemicellulose A, the molecular ratio of these compounds is 10:1, while 
in B, it is 5:2. Under conditions of partial hydrolysis, both the fractia 
give rise to glucose and a gluco-galacturonic acid. 
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